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THERMAL-HYDRAULICS MODULE, TH-2 

LIQUID METAL FAST BREEDER REACTOR THERMAL-HYDRAULICS 

, ^ *" - ■* ' 

1.0 Oblect^f 'Module - 

The objeibt of this^module is to present the basic e •u'r f liquid 

metal fast breeder reactor thermal-hydraulics. This reo' ' 
demons tfa1:ion of: 



• How the actual fuel geometry can be modeled for simplified 
thermal-hydraulic analysis. 

• What information is necessary to characterize the t 'armal- 
hydraulic behavior of the reactor. 

• The development of the theoretical relations that permit 
the computation of these thermal-hydraulic characteristics. 

• The actual calculation of this information for the reactor. 
This calculation requires the use of the LMFBR Thermal- 
Hydraulics Code, the description of which is included in 
this module. 

The thermal-hydraulic characteristics of the reactor are required for 
the determination of: 

• Fuel integrity 

• Cladding integrity 

• Coolant exit conditions 

• Pumping requirements 

• Temperature feedback for reactor neutronics calculations. 
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2 . 0 Content of the Module 

This learning module contains the thermal-hydraulics of liquid metal 
fast breeder reactors. Specifically, the module is concerned with the 
temperature field, the heat transfer rates and the coolant pressure drop 
in typical LMFBR fuel assemblies. 

As in all of the modules of this series, emphasis is placed on developing 
the theory and demonstrating its use with a simplified model. The model is 
carefully selected to insure that analyses based on it will exhibit all of 
the important thermal-hydraulic trends of the typical reactor. The 
description of the fuel for a typical LMFBR and the modeling of the fuel 
which preserves the significant thermal-hydraulic characteristics are treated 
in the next section of this module. 

Following the geometry and modeling section, the basic theory governing 
the temperature distributions, heat transfer rates, pressure drops, and 
energy balance considerations is presented. The temperature distributions 
in the fuel and cladding are calculated assuming one-dimensional radial heat 
conduction. The pressure drop in the coolant channels and the heat transfer 
coefficient for use in Newton's law of cooling are calculated from empirical 
relations developed for reactor coolant channel flows. Energy balances for 
small axial segments of the coolant channel are used to step the solution in 
the axial direction. Simple examples, illustrating the individual 
calculations, are worked out in detail for typical LMFBR conditions. 

The heart of the module is the LMFBR Thermal-Hydraulics Computer Code. 
Basically, the code solves for the radial temperature distributions in the 
fuel, cladding, and coolant at any axial station and then marches axially 
with an energy balance in the coolant. The code and its use are described 
in detail. Included are a listing and definition of all variables, a 
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discussion of all input require:nents and resulting output, an annotated 
flow chart of the code, an explanation of all options in the code, and a 
listing of the code which includes enough comment statements to clearly 
indicate the operational steps being performed. By proper specification of 
the options, the code can either be used as an individual entity to study 
thermal-hydraulic aspects exclusively, or as a subroutine in the total LMFBR 
module package to provide temperature feedback to the other modules. 
Examples are worked out using the code. In typical examples, the location 
and magnitude of the maximum fuel temperature in the LMFBR are found and the 
effect of oversized cladding on maximum fuel temperature and coolant outlet 
temperature are determined. 
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3.0 LMFBR Fuel Geometry and Its Model 

The typical LMFBR core is approximately cylindrical in shape and is 
composed of about 200 to 250 fuel assemblies. These fuel assemblies are 
hexagonal in cross section, 10 to 15 feet long, and stand vertically in 
the core. Small gaps separate the fuel asseniblies from one another. The 
fuel assemblies are contained inside a core shroud which in turn fits inside 
a cylindrical core support which is 10 to 15 feet in diameter. To provide 
a flattening of the neutron flux (and hence the power production) in the 
radial direction, the enrichment of the fuel is varied in cylindrical zones. 
The individual fuel assemblies are orificed such that the coolant flow is 
proportional to the power generation of the assembly. 

The fuel assenibly consists of 150 to 275 fuel rods, support grids 
holding the rods firmly at the ends, and a hexagonal flow duct. The spacing 
between the individual rods and between the rods and the duct wall is main- 
tained by helically wound wire wraps on the fuel rods. About 4 to 6 feet of 
the fuel rods contain nuclear fuel and therefore generate heat internally 
during reactor operation. 

A cross-section of a typical fuel assenibly is shown in Figure 1. The 
fuel rods are arranged in a triangular array within the fuel assembly. The 
liquid metal coolant, which carries away the energy generated within the fuel 
rods, flows in the spaces between the rods and between the rods and the duct. 
Three types of coolant channels can be identified as shown in Figure 1. These 
correspond to the channels associated with central fuel rods, outer row fuel 
rods, and corner fuel rods. 

Note that a great majority (384 out of 438 coolant channels in a 217- 
pin bundle) of the coolant channels are central-type coolant channels. Also 
note that there are about two (actually 2.02) flow channels for each fuel rod. 

ft 
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Figure 1. LMFBR Fuel Geometry. 
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Thus, to a reasonable approximation, the coolant passage associated with 
each fuel rod in the model can be taken as two central-type channels. The 
model flow passage per fuel rod, therefore, consists of two triangles, 
each with corners at the centers of the three fuel rods forming the channel 
minus the three sixty-degree segments of the fuel rods themselves. The 
ratio of center-to-center spacing, S, to the fuel rod diameter, D^, is 
typically 1.1 to 1.3. 

Each flow channel serves segments of three fuel rods as can be seen in 
Figure 1. Since the three rods are identical, and since the neutron flux 
does not vary significantly in the short distance from one rod to the next, 
all three rods may be assumed to have the same operating characteristics. 
Thus the three fuel rod segments associated with each flow channel can be 
treated as one-half equivalent fuel rod. The power generation level in this 
model fuel rod will be taken as that of the average fuel rod in the core. 

The model fuel rod will have a variation in power density in the axial 
direction. This variation will be taken as 

q (Z) = q cos 

e 

unless a computed actual axial variation is provided from another module. 
In Equation 1, q' " represents the thermal source strength per unit volume 
at any axial location Z, q'^' represents the thermal source strength per 
unit volume at the center of the fuel rod (Z=0) , and is the extrapolated 
height of the core. Both q"' and q'^' are taken to be constant radially 
throughout the fuel in the fuel rod. The magnitude of q''' is representative 
of that of the average fuel rod in the core. The extrapolated height is 
calculated from 

H = H + 2 L (2) 

' id 
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where H is the actual height of the reactor core and is an extrapolation 
length. The extrapolation length is the distance between the actual end 
of the core and the location where an extrapolation of the waveform repre- 
senting the actual neutron flux distribution within the core goes to zero. 
From neutron diffusion theory can be shown to be about one migration 
length for a bare core. The migration length can be calculated from 
neutron diffusion theory [1,2,3], 

The LMFBR fuel rods are 10 to 15 feet long and about 0.25 to 0.30 inch 
in diameter. About 4 inches of the total length is used for end fittings, 
about 1 foot on each side of the active core consists of a blanket of depleted 
UO2, and about 3 to 5 feet near the downstream end serves as a fission gas 
plenum. The remaining length contains nuclear fuel, usually uranium dioxide 
and Plutonium dioxide. The fuel is typically in the form of cylindrical 
pellets with a diameter of about 0.25 inch. These fuel pellets are stacked 
inside a tube of cladding, a stainless steel alloy with about a 0.02 inch 
thick wall. There is a 0.003 to 0.005 inch gap between the fuel and the 
inside of the cladding. 

For simplicity, the gap is neglected in the model of the fuel rod. 
The model fuel rod thus consists of fuel of diameter contained in a 
tube of cladding with an inside diameter and an outside diameter D^. For 
the model fuel rod, the length of the end fittings will each be assumed to 
be 3% of the active fuel length, the length of the blankets will each be 
assumed to be 20% of the active fuel length and the length of the gas plenum 
on the top of the fuel will be taken as 80% of the active fuel length. 



^Numbers in brackets refer to items in References. 
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Each fuel rod is wound with a wire wrap to provide spacing and also 



to promote mixing between the coolant channels. For the typical S/D^ 



ratio of 1.2, the diameter of the wire wrap corresponds to 20% of the 
outside diameter of the cladding- The helical pitch with which the wire 
is wound on the rod is typically 1 foot. 

The wire wrap on the fuel pins complicates the calculation of the 
equivalent diameter of the coolant channels. The wire wrap both reduces the 
flow area and adds to the wetted perimeter. In making the area and wetted 
perimeter calculations, it is important to note that there is a wire wrap 
inside a given flow channel for exactly 1/2 of the total axial distance. 
Thus it is reasonable to treat each flow channel as if 1/2 wire wrap were 
present at each axial level. The flow area and wetted perimeter then are 



From these, the equivalent diameter of the coolant channel is given as 



= y cos 30 




(3) 



p = ^ (D + D ) . 
*^w 2 s WW • 



D = 




(5) 
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4.0 LMFBR Thermal-Hydraulic Theory 

4.1 Internal Heat Generation 

As a result of nuclear fission in the fuel, heat is generated. The 
rate of energy generation in the fuel per unit volume is called the "volu- 
metric thermal source strength*', q''', and can be calculated from 

f 1 1 — 

" ^f^ff ""f ♦ (6) 

where is the energy released per fission reaction (energy dimensions), 
^f f fissionable fuel density (fissionable nuclei per unit volume) , 

is the effective fission microscopic cross section (dimensions of area), and 
(j) is the neutron flux (neutrons per unit area per unit time). Note that q''' 
has dimensions of energy per unit volume. The details of this calculation 
of volumetric thermal source strength are found in Reactor Statics Module 8. 

The volumetric thermal source strength varies throughout the reactor 
since (J), and perhaps N^^, vary. For a cylindrical reactor, axial symmetry 
is generally a reasonable approximation. Therefore q''' reduces to a function 
of only the radius, R, and the axial position, Z. Across any single fuel 
rod, the change in q'*' is small (due to small change in R) and q''' can 
therefore be considered constant across its cross section. Thus for any 
single fuel rod q'*' is a function of only the axial position. Also for any 
small axial segment of a fuel rod (the order of 1% of the total active length) 
the change in q*'' is moderate. Using a constant average value for any 
such segment thus introduce little error into the analysis. This approximation 
is used in this module. 

The heat generation in any small segment of the fuel can be obtained by 
multiplying the volumetric thermal source strength by the volume of the 

o 13 
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segment. For steady state conditions this energy must be removed from 
the fuel. The mechanism by which this heat transfer occurs within the 
fuel is thermal conduction. 

Example 1. 

A LMFBR containing 25,000 fuel rods generates 2000 Mw (t). The 
cylindrical core of the reactor is 5 ft long and 4 ft in diameter. The 
fuel consists of 0.22 in. diameter fuel pellets. These are contained 
in cladding which has an outside diameter of 0.25 in. The fuel rods are 
spaced 0.30 in. center- to-center in a triangular array. The volumetric 
thermal source strength of the fuel varies as q'" = q'^' cos(nZ/H^) . 
The extrapolation length for the fuel is 3 in. Find the power per unit 
length (in kw/ft) and q'^' for the average fuel rod in the core. 
Solution . 



The power generated in the average fuel rod is 

^ Ootal ^ 2000 X 10^ „ 80.000 watt = 80 kw. 
ave N 25,000 

. , , P 80 kw kw 

Thus, per unit length the power is jj = ~3~f7 = 1° J^- 

The power produced by the average fuel rod can also be calculated 
by Integrating the power produced in each differential volxme of the fuel 
rod. Thus, 

fH/2 



P 

ave 



q' "dV 



q'V cos (nZ/Hj(nD^/4)dZ 

i-E/2 



(H q'''/4)[sin(nZ/H )] 
e r o e 



= (H^ q'o'/2) sln(IIH/2H^)u 
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Solving for q'^': 



2P 

. axe 



H d" sin(nH/2H ) 
e f e 



(2) (80,000 watt) (3.413 ^^^t-hr^ 



(5.5 ft)(0.22 in.)^(-j2f^.>^ sin[(n)(5 ft)/(2)(5.5 ft)] 



2.98 X 10^ Btu/hr-ft"^. 
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4.2 Temperature Distribution In the Fuel 

For conduction heat transfer the heat flux Is proportional to the normal 
temperature gradient. When the proportionality constant Is Inserted, 

ST 

where is the heat transfer rate in the n-direction (energy/time) , 8T/8X^ 
is the temperature gradient in the n-dlrectlon (degrees/ length) , is the 
area normal to the n-direction, and k is the thermal conductivity of the 
material (energy/degree-length-time). Equation 7, which is Fourier's 
law of heat conduction, relates the heat transfer to the temperature field 
and is also the defining equation for the thermal conductivity. The thermal 
conductivity is a material property and its magnitude in general varies with 
the temperature of the material. Heat transfer properties of various 
reactor materials are tabulated as functions of temperatures in References 
4 and 5. To reduce the complexity of heat transfer calculations, the thermal 
conductivity is often assumed to be constant and evaluated at an average 
temperature. The minus sign in Equation 7 assures that the heat transfer 
is in the direction of decreasing temperature. Equation 7 shows that 
temperature gradients are required for heat transfer. In nuclear reactor 
applications, where there are high heat transfer rates, large temperature 
variations occur. One of the primary tasks of reactor thermal-hydraulic 
analysis is the prediction of this temperature distribution in the fuel. 
The heat conduction equation in cylindrical coordinates, 

.32t, 1 3T . l_3!T^i!T^£lll. 1 |I 
-;^^rl^^-^—2\^2^ ^ a . (8) 

and the initial and boundary conditions prescribe the temperature distribution 
within the fuel rod. In this relation, the thermal conductivity, k, and the 

i6 
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thermal dif f usivity , a, have been assumed constant. The development of this 
equation is found in heat transfer texts [4,5,6] and several simplified 
cases are given as exercises for the student. For steady state conditions, 
the unsteady term on the right hand side of the equation is zero, and for 
axial symmetry, the 0 variation disappears. It may also be observed that 
since the length of a fuel pin is much greater than its radius, the tempera- 
ture gradients in the radial direction will be much greater than the tempera- 
ture gradients in the axial direction. Therefore, to a good approximation, 
the heat transfer in the axial direction can be neglected with respect to 
that in the radial direction, and the resulting temperature distribution 
and heat transfer reduced to a one-dimensional case for any axial segment 
in which q''' may be assumed constant. Th*^ governing differential equation 
for this case reduces to the ordinary differential equation 

d^T^ldX + alll -.0; (9) 

The solution of Equation 9 subject to the boundary conditions 

T = T at r = 0, 
o * 

dT (10) 

-7^ - 0 at r = 0 

dr 

yields the temperature distribution in the fuel. The second boundary 
condition is obtained from the observation that the temperature distribution 
must be continuous across the center of the cylinder. Observing that both 
q'*' and k are constant and that 

dr2 

Equation 9 can be written as 



,0 r dr " r dr I dr / 



(11) 
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Integrating twice results in 

rp . . li + Cj^ In r + C^. (13) 

Applying the boundary conditions of Equation 10, the integration constants 



are 



= 0, 
S = ^o- 



(14) 



Substituting into Equation 13 gives 



.tit 



T - To " 4k 



f 



(15) 



This relation shows the teinperc>f;ure distribution in the fuel to be parabolic 
with inaxlnnim temperature at the center. The heat transfer rate through any 
cylindrical shell can be calculated from Fourier's law which takes the form 



<»r - - '^f ^ f (^^^ 



where A = 2Trr(AL) , AL being the length of the cylindrical shell, 
r 

Of particular interest are the temperature and heat transfer at the 
surface of the fuel; i.e., at r = = D^/2. At this location. 



q'*' 2 

■ " 4kf • ""f » (17) 



(18) 

.1 t I 



dTI 

- -kj2irr^ (AL) a^rlr - r, 

- ifrj2(AL) q' 

Noting that Trr^AL is the volume of the fuel rod segment of length AL it is 
observed that the heat transfer out of the r = r^ cylindrical shell is indeed 
equal to the total energy generated as calculated from (q'") (fuel volume). 
For one-dimensional heat transfer all of the energy generated within the 
fuel must be transferred out through the surface. 

18 
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Example 2.. ■ ^j. 

Equation 17 relates the temperatures at the center and surface of a 

a * - 

i * . . ^ 

heat-generating cylindrical fuel rod. This equation contains the thermal 
conductivity of the fuel which, in general, is a function of temperature. 
However, the analysis leading to Equation 17 assumes the thermal conduc- 
tivity to be constant. Determine the temperature at which to evaluate 

to make the assumption of constant k^ compatible with a k^ that varies 
linearly with temperature. 
Solution. 

The method of solution is to develop an expression relating T^ and 
T^ which assumes k^*= a + bT and then compare this result with Equation 17. 
We begin with Fourier's law of heat conduction, 

^ = -^f^rd?- 



Recognizing that q = q''' (volume) and = 2iTr(AL); 



q"' TTr^(AL) = k^ 2Trr(AL)-^. 



Simplifying, separating variables, and substituting for k^: 

«• • X' 

q'" rdr = -2(a + bT) dT . 
Integrating from r = 0 to r = r^ with q'*', a, b constant, 



= -2 



0 

2 

■J. _ 



aT + lr^ 
l_ _lT 



-a(T^ - T„) - |(T'f - T^)2 



= [a+f(T^ + Tp](T^-Tp 
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T - = 
o f 



Comparing this with Equation 17 written as 

III 2 i 

q 

o f 4 k- 



we see that the relations are equivalent for 

a + |(T^ + Xp. 

Note that this is precisely equal to evaluated at the arithmetic average 

T + 

I/O tx 

fuel temperature; i.e., k( 2 
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4.3 Temperature Distribution In the Cladding 

The temperature distribution in the claddlnR can be calculated by 
methods similar to those used for the fuel. Noting that uo energy is 
generated in the fuel, the governing equation simplifies to 

|j (r f) . 0. an 

The boundary conditions obtained by matching temperature and heat transfer 

at the fuel-cladding interface are 

T « at r ■ r^, 

dT ^f ^f^'" (20) 

dr " " 2TTk r.(AL) " " 2k at r - 
c I c 

where k is the thermal conductivity of the cladding, 
c 

Integrating Equation 19 twice yields 

T - C3 In r + C^. (21) 

Applying the boundary conditions in Equation 20, 

S ■ " ~2ir ' (22) 

Substituting into Equation 21 yields the temperature distribution in the 
cladding. 



f ^ 1 r 



Or, if q''' is eliminated in favor of q^. 



(23) 



The temperature at the surface of the cladding, T^, is obtained by evaluating 

Equation 24 at r , 
s 
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T = - ^ , In — . (25) 

s f 2iTk (AL) r-: 

c t 



Note also that the same heat transfer rate must occur through each 
cylindrical layer for the steady state conditions ; i . e . , = q^- 



S2 
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4.4 Temperature Drop In the Convectlve Layer 

All of the heat generated within the fuel ultimately must be transferred 
to the coolant. Thus, all of this energy mist be transferred through the 
fluid layer near the outer surface of the cladding. The heat transfer 
mechanism, wherein the energy is carried away or convected from the solid 
surface by a fluid in motion, is called convection. The heat transfer rate 
for convective heat transfer is related to the temperature difference 
between the surface and the bulk fluid, the driving potential for the heat 
transfer, by Newton's law of cooling, 



This relation may also be taken as the defining equation for the heat transfer 
coefficient, h. The bulk temperature, T , is a mass-weighted average 

D 

temperature of the fluid in the flow channel. It is formally defined by 



This is the temperature that a thermometer would indicate if immersed in the 
mixed fluid collected from the discharge of the flow channel at the given 
location. 

Relations for calculating h from the flow characteristics and coolant 
properties have been developed. These will be treated in a later section of 
the module. For the moment, the quantity of interest is the temperature 
drop across the convective layer. In terms of h this can be obtained from 
Equation 26 as 



(26) 




T„ - T. 



2ifr (AL)h 



(28) 



where r is the outer radius of the cladding. 



s 
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Example 3> 

Consider the LMFBR described in Example 1. The bulk temperature of 
the coolant changes from 665 F to 1335 F as the coolant flows through 

the reactor. Given that the heat transfer coefficient is 44,000 

2 

Btu/hr-ft -F, the thermal conductivity of the cladding is 12.45 
Btu/hr-ft-F, and the thermal conductivity of the fuel is 2.00 Btu/hr-ft-F, 
determine the fuel centerline temperature at the axial center of the core. 
Solution. 

Half of the total energy transfer to the coolant occurs in the first 
half of the core. Thus the coolant bulk temperature at the axial center 
of the core is 1000 F. To this, we must add the temperature rises across 
the convective layer, the cladding, and the fuel to obtain the fuel 
centerline temperature. 

The heat transfer rate through the convective layer at Z « 0 is 

q = q'" (fuel volume) = q'*' ttD^(AL)/4 
o ^ o r 

= 2.98 X 10^ Btu/hr-ft"^)(7T/4)(0.22 in.)^(ft/12 in.)^(AL ft) 
= 78,784 (AL) Btu/hr . 

From Equation 26 , 

m 

^ ^ % 78784 AL Btu/hr 

J _ ^ 5= - 

s ^ h (44,000 Btu/hr-ft -F)7r(0.25 in.)(ft/12 in.)(AL ft) 

» 27.3 F. 

The heat transfer rate through the cladding is the same as that 
through the convective layer. Equation 25 yields 
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' " 2Tik (AL) r, 
c f 



_ (78784 AL Btu/hr) In (0. 125/0. 11) 
2Tr(12.45 Btu/hr-ft-F)(AL ft) 



= 128.7 F . 

From Equation 17 the temperature rise through the fuel is found as 



1 ^ 2.98 X 10"(Btu/hr-ftn(Q.ll in.) (ft/12 in.)' 



~ ^f 4 4(2.00 Btu/hr-ft-F) 

= 3130.0 F . 

Adding these temperature rises to the bulk temperature, 



= (T„ - + (T^ - T^) + (T^ - T3) + 



= 3130.0 + 128.7 +27.3 + 1000 = 4280.0 F. 
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4.5 Pressure Drop In the Coolant 

All of the energy generated by nuclear fission in the fuel must be 
carried out of the reactor by the coolant which flows in the channels 
between the fuel rods. In LMFBR the coolant is a liquid metal which is 
force-circulated by a pump. Enough pressure head must be provided by the 
pump to overcome the pressure losses in the flow loop. The entire loop 
consists of the reactor flow channels, the heat exchangers, coolant purity 
maintenance components and connecting piping. Of interest here is the 
determination of the pressure loss incurred in the flow through the reactor 
core. Two types of pressure losses will be considered. These are the 
frietional pressure loss along the coolant channel and the entrance and 
exit permanent pressure losses. 

The frictional pressure loss is a manifestation of the shear stress 
on the flowing fluid by the walls of the flow channel. This pressure loss 
is calculated from 



This relation may also be interpreted as the defining equation for the 
friction factor, f. The bulk velocity, U^, found in Equation 29 is the 
velocity averaged across the flow channel. It is defined by 



AP 




(29) 




'A 



(30) 



This velocity is related to the mass flow rate by 



m = pAU. 



(31) 



B 



or to the mass velocity by 



G =^= pUg. 




(32) 



The equivalent diameter, D^, is used to relate experimental data from 
circular tubes, which is abundant in the literature, to different geometries. 
The equivalent diameter is calculated from 

D = . (5) 

2 

Note that for a circular cross section the flow area is ttD /4 and p = ttD, 

'^w 

Thus for a circle the equivalent diameter is equal to the actual diameter, 
as required. This, of course, is why the definition of Equation 5 is used. 

The value of the friction factor depends on the flow and surface 
conditions in the channel. The flow conditions are characterized by the 
Reynolds Number, 

PU-D 

Re = (33) 

y 

It is well known that for Re below a critical value the flow is laminar, 
and for Re above the critical value the flow is turbulent. For internal 
flows this critical Re is about 2000. Laminar flow may be thought of as an 
ordered process in which fluid layers slide over one another, being retarded 
only by the molecular interaction between the layers. The viscosity of the 
fluid quantifies the magnitude of this interaction. In laminar flow any 
disturbance in the fluid is damped by the viscous action. In turbulent 
flow there is an additional random transport mechanism operable. This 
mechanism may be modeled as eddies (finite sized patches of fluid which 
retain their characteristics for finite times) moving throughout the fluid, 
transporting mass, momentum, and energy by virtue of their movement. This 
action is quite violent and results in transport rates much greater than 
those by purely molecular activity in laminar flow. In turbulent flow, 
disturbances in the flow field grow and propagate resulting in increased 



turbulence levels downstream of the disturbances- Even in highly turbulent 

Clows, there exists a layer near the solid bounding surfaces where the 

presence of the wall retards the penetration of eddies and therefore acts 

as a flow stabilizer. This results in a laminar sublayer existing near any 

solid boundary. Even though this layer is very thin, it is very important, 

since much of the temperature and velocity change between wall and bulk 

conditions occurs in this layer. In reactor coolant channels the flow is 

highly turbulent and Re of the order of 100,000 arc quite common. i ^- 

The surface conditions of the flow channel is characterized by the ratio 

of the surface roughness height to the equivalent diameter of the flow 

channel, e/D . The flow channel exhibits smooth tube behavior if the rough- 
e 

ness height is less than the thickness of the laminar sublayer. For reactor 
coolant flow channels the surface conditions are controlled so that this 
criterion is met. 

The calculation of the friction factor for the reactor coolant flow 
channels is a three-step procedure. First the smooth circular tube friction 
factor is determined, then a correction factor is applied to account for 
the actual channel shape, and then a correction factor is applied to account 
for the wire-wrap spacers. The smooth circular tube friction factor can be 
determined from a Moody chart where f^^.^.^ is plotted versus Re. These charts 
are found in most fluid mechanics texts and handbooks [7,8,9,10]. An alter- 
nate method is to calculate the friction factor from a correlation in equation 
form. One of the most widely accepted for reactor work is 



f = 0-184 . (34) 
circ - ^^0.2 • 



This relation, like the Moody diagram, was obtained from a curve fit of 
experimental data in smooth circular tubes. The friction factor defined in 
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this module is the Darcy-Weisbach friction factor. Care must be taken 

not to confuse this friction factor with the Fanning friction factor (^Yan " 

Tt/^) also found in the literature. 
D-W 

A correction of the circular tube data to the rod bundle geometry has 
been obtained by Deissler and Taylor [11]. Their results, expressing ^o^^circ 
as a function of S/D^, can be approximated by the following polynomial fits 



f /f . - -3.475 + 8.053 |- - 4.705 (|-)^ + 0.9162 (35) 

s s s 



The correction for the wire-wrap spacer on the friction factor is quite 
significant. This correction may increase the smooth channel value by the 
order of 30% for typical LI^BR applications The most extensive pressure 
drop data on wire-wrapped reactor fuel bundles in LMFBR operating ranges 
is that of Reihman [12]. His data show the effects of wire-wrap size, 
channel equivalent diameter, and wire-wrap helical pitch. The results may 
be approximated in polynomial form by 

D D 2 

f/f . 1.195-(HP-0. 75)10.554-5. 456 :ff5^+ 11. 58(^) ]. (36) 
o e e 

The permanent pressure losses at the inlet and exit of the individual 
flow channels are the result of increased viscous energy dissipation resulting 
from increased turbulent activity. These losses may be calculated from 

»E ■ K ^ . (37) 

In this relation K is the resistance coefficient for the inlet or exit 
geometry. The inlet may be approximated as a sudden contraction from a 
very large diameter to the equivalent diameter of the flow channel. For this 
sudden contraction, the resistance coefficient is 0.5 [10]. Similarly the 
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exit may be approximated as a sudden expansion from to a very large 
diameter. For this expansion, K is 1-0 [10]. 

For each fuel assembly there is an additional pressure loss due to 
orif icing. This additive pressure loss is omitted in this module, not 
because it is unimportant, but because its calculation requires specific 
design information on the individual fuel assemblies. 
Example 4. 

The LMFBR described in Example 1 has an inlet coolant temperature of 
665 F, an inlet coolant pressure of 200 psia, and a coolant mass velocity 
of 7.0 X 10 Ibm/hr-ft . The coolant exit temperature is 1335 F. 
Estimate the frictional pressure loss in the core. 
Solution. 

The coolant properties are evaluated at the average coolant tempera- 
ture. At 1000 F: 

VI = 0.505 Ibm/hr-ft, 
p = 51.19 Ibm/ft^ . 
The average coolant velocity is found from 

U = G/p = (7-0 X 10^ lbm/hr-ft^/(51-19 Ibm/ft^) = 136,800 ft/hr . 
B 

The equivalent diameter is computed next. 

Af = ^ S2 cos 30° - f(Dj + D^) 

= 1/2(0 . 300 in.)^(0.866) - |[ (0.250 in.)^ + (0.050 In.)^] 

. 2 
= 0.0134 xn. 
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p = ^(D + D ) 
•^w 2 s WW 



= |-[ (0.250 in.) + (0.050 in.)] = 0.471 in, 



D = 4 A-/p = 4(0.0134 in.^)/(0.471 in.) 
e r w 



= 0.114 in. = 0.00948 ft 



The Reynolds niunber is 
Re = Ug p/v. 

= (136 , 800 ft/hr) (0.00948 ft) (51. 19 Ibm/f t'') /(0 . 505 Ibm/hr-ft) 
= 131,400 . 

Using a circular tube friction factor correlation, 

f , = 0.184/Re°'^ = 0.184/(131,400)°'^ 
circ 

= 0.0174 . 

Applying the correction for the actual reactor flow channel, 

f = f . [- 3.475 + 8.053(S/D ) - 4.705(S/D )^ + 0.916(S/D )^] 
o circ s s s 

= (0.0174) [- 3.475 + 8.053(^^) - -^.705(1^)^ + 0.916(^y||)^] 

= 0.0183 . 

Applying the correction for the wire wrap , 

f = f [1.195 - (0.186 - 1.364(D /D ) + 2.894(D /D )^)(HP- 0.75 ft)/0.25 ft] 
a • WW 3 WW e 
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- (0.0183)[1.195 - (0.186 - 1. 364(§^551^) + ^.ml^^^^) ) 

(1 ft - 0.75 ft)/(0.25 ft)] 

« 0.0192 • 

From the defining equation for the friction factor, 
APp - f(AL/D^) p Ug/2g^ 

= (0.0192) (60 in. /0. 114 in.) (51. 19 Ibm/ft^) (136,800 ft/hr)^ 
/(3600 sec/hr)^(2)(32.2 Ibm-f t/sec^-lbf) 

= 11600 Ibf/ft^ 

= 80.5 psi . 
Example 5. 

Each coolant channel in a fuel assenibly has the same pressure drop. 
If all coolant channels are not identical the flow will redistribute such 
that this equal pressure drop is attained. Determine the relation between 
average velocity in the coolant channel (or equivalently the mass velocity) 
and the equivalent diameter of the channel which governs the flow distri- 
bution. Neglect the effect of the equivalent diameter in the wire wrap 
friction multiplier. 
Solution. 

The frictional pressure drop is given by 

e c 
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Substituting for the friction factor, 



0.184 AL P"b 

Re°-2D 2g 
e c 



0.184 AL P^^B 

„ „ 0.2 D 2g 

U„D 0 e c 
B e 



D e 



if constant density and viscosity are assumed. For channels 1 and 2, each 
having equal pressure drop , 



or 



"b2/"b1 = (°e2/°el>'^' 
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4.6 The Heat Transfer Coefficient 

The usual relations for heat transfer coefficient developed for water, 
air, and similar "normal" substances can not be used for liquid metals. 
The reason is that the Prandtl number range of liquid metals (typically 0.001 
to 0.01) is vastly different than that of the other substances (typically 
0.7 to 30). Consequently, analyses based on the assumption that Pr = 1 and 
correlations obtained from data in the 0.7 to 30 range of Pr could not be 
expected to yield even order of magnitude results. 

The Prandtl number is defined as the ratio of the molecular diffusivity 
of momentum to the molecular diffusivity of heat; i.e., Pr = v/a where v = 
y/p and a = k/pC^. The Pr of liquid metals is very low compared to that of 
other fluids because of the very high thermal conductivities (typically 40 
Btu/hr-ft-F compared to 0.1 Btu/hr-f t-F) . Tne Prandtl number relates 
the temperature field to the velocity field. The low value of Pr for liquid 
metals indicates that the velocity field is relatively unimportant in the 
establishment of the temperature field in the fluid. The low value of Pr 
also indicates that the resistance to heat transfer is now distributed more 
uniformly through the fluid rather than largely concentrated in the laminar- 
sublayer as it is for higher Pr fluids. Thus, it should be apparent that 
liquid metal heat transfer is significantly different than that in ordinary 
fluids . 

For turbulent flow in circular tubes with uniform heat input several 
correlations appear in the literature [13,14,15,16]. Of these, perhaps the 
best known is the Lyon-Martinelli relation [13,16], 

Nu - 7 + 0.025 Pe°-' • 
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The two dimensionless moduli in the equation are the Nusselt nuniber, 




(39) 



and the Peclet number, 



Pe = Re Pr 



(40) 



k ""a 



One common interpretation of the Nu is that it is the ratio of the actual 
heat transfer from a surface to the heat transfer that would result by mole- 
cular conduction through the stagnant fluid. The Pe, rather than the Re 
and Pr individually, is common to all liquid metal heat transfer coefficient 
correlation equations. This grouping eliminates the viscosity from the 
correlation and thereby signifies that the viscosity is not rate-determining 
in liquid metal heat transfer. The theoretical basis of the Lyon-Martinelli 
relation stems from the analogy of heat and momentum transfer by Martinelli 
[16]. Lyon [13] assumed that the turbulent eddy conductivity of heat is 
negligible compared to the high molecular conductivity of liquid metals 
and obtained the simplified result of Equation 38. 

Another common liquid metal heat transfer correlation for circular tubes 
is that of Lubarsky and Kaufman [14], 



This relation was obtained by carefully analyzing and reevaluating all the 
existing liquid metal heat transfer data at the time an*i then obtaining the 
best fit to the data. This relation predicts considerably lower values of 
Nu than Equation 38. Since this is the more conservative result, the use 
of Equation 41 rather than Equation 38 is recommended. It should be 
pointed out however, that the prediction of he?t transfer coefficient for 



Nu - 0.625 Pe 



0.4 



(41) 
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liquid metal systems is not a precise science • Even after rejecting some 
data for oxide contamination, the remaining data shows deviations of up 
to ±50% from the values predicted by the correlations. 

Heat transfer coefficient data for flow parallel to rod bundles is 
limited- The only experimental data in the LMFBR range of S/D^ is that of 
Borishansky and Firsova [17] and Subbotin et al. [18] and there are serious 
questions as to oxide contamination and surface wetting problems in these 
data. Also, since the velocity field and temperature field are not strongly 
coupled for low Pr fluids, the multiplicative friction factor correction of 
Deissler and Taylor [11] does not apply. Therefore reliance will be placed 
on the analytical work of Dwyer [19] in which he predicts Nu for liquid metals 
in fully developed turbulent flow in staggered rod bundles with uniform 
heat flux. Fortunately, Dwyer 's results are directly applicable to the 

low S/D range (1.1 to 1.3) of the LMFBR. There is some experimental verifi- 

s 

cation of Dwyer 's predictions in that they agree well with the data of 
Borishansky and Firsova for S/D^ =1.2 and Pe greater than 200. The final 
results of Dwyer 's work are somewhat limited due to the lengthy computational 
time required. The results that are available can be approximated in the 
following form which is convenient for use in the LllFBR Thermal Hydraulics 
Code : 



Nu 



for 



Pe i Pe 



o 



o 



Nu = 



0.8 



(42) 



NUq + 0.025(Pe-PeQ) 



for Pe > PCq^ 



where 




(43) 



Pe = 336 |- - 110 . 



(44) . 



36 



- 33 - 



This approximation, in which the equivalent diameter is to be used in the 
Nu and Pe, is expected to fit the results of Dvjyer within 15% in the range 
of S/D^ from 1,1 to 1,5 and the Pe range from 30 to 3000. 

The above results apply strictly only to staggered rod bundles without 
wire wraps. There is no data in the literature on the heat transfer 
characteristics of wire-wrapped bundles • One could postulate various effects 
of the wire wraps on the heat transfer coefficient • For instance, the 
increased turbulence from the wire wraps would tend to increase the heat 
transfer. But for liquid metals where molecular conduction far outweighs 
eddy conduction one would expect this effect to be secondary. The wire wraps 
also add appreciable heat transfer area (S/D^ times the unwrapped area). 
However, any heat transfer into the liquid metal through the wire wrap would 
experience resistance to heat transfer in the wire wrap itself and a contact 
resistance between the wire wrap and the cladding. Finally, there is an 
effect via the equivalent diameter, the wire-vrapped bundle having a sub- 
stantially smaller equivalent diameter. Thus, when the heat transfer 
coefficient is calculated from the Nusselt number, a larger value of h will 
result. There is also a secondary effect in the opposite direction in that 
Pe will be smaller. The effect just described via the equivalent diameter 
can be justified physically. The smaller flow channel has a lower thermal 
resistance between the wall and an "average" fluid particle in the channel 
(since the mean distance is smaller) . The procedure recommended and the one 
used in the LMFBR Thermal-Hydraulics Code is to use the equivalent diameter 
calculated with the wire wrap in place and neglect the other two possible 
"corrections" on the heat transfer. 
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Example 6> 

Determine the average heat transfer coefficient for the LMFBR 
described in Examples 1 and 4, 
Solution, 

The average heat transfer coefficient corresponds closely to 
that calculated at the average coolant conditions; i.e., 1000 F and 
160 psia. Additional fluid properties required for this calculation are 

C « 0.312 Btu/lbm-F, 
P 

k = 37.62 Btu/hr-ft-F . 

The Prandtl number is 

Pr « C y/k 
P 

= (0.312 Btu/lbm-F) (0.505 Ibm/hr-f t) /37.62 Btu/hr-ft-F) 

= 0.00419 . 

Computing the Peclet number , 

Pe = Re Pr « (131,400) (0.00419) = 550.3 - 

Applying the Nu correlation for reactor flow channels > 

Pe = 336.0(S/D ) - 110.0 
o s 

= 336.0(0.300 in. /O. 250 in.) - 110.0 
= 293.2 

Nu = - 88.0 + 130.8(S/D ) - 4l.67(S/D )^ 

O 8 S 
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= - 88.0 + 130.8(0.300 in. /O. 250 in.) - 41.67(0.300 in. /O. 250 in.)^ 
= 8.955 

Nu = Nu + 0.025(Pe - Pe )°'^ 

o o 

= 8.955 + 0.025(550.3 - 293.2)°'^ 
= 11.073 . 

Calculation of h from the definition of Nu gives 
h = (Nu) k/D 

e 

= (11.073) (37.62 Btu/hr-ft-F)/(0.009'+8 ft) 
= 43,940 Btu/hr-ft'^-F • 
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4.7 Incremental Energy Balances 

As heat is transferred from the fuel rod to the coolant, the temperature 

of the coolant rises. Thus, the coolant temperature varies from a minimvsm 

at the inlet of the coolant channel to a maximum at the exit of the fuel 

bundle. The bulk temperature at any axial position can be determined from 

an energy balance where the total energy added by heat transfer to the coolant 

is equated to the energy rise of the coolant; i.e., 

Z 

ii - il , + / (Z) dZ. (45) 

xn J • 

-H/2 

Recognizing that for liquid metals Ai - AT^ and solving for the local bulk 
temperature, this becomes 



Equation 46 can be used to determine the bulk temperature Increase for 
any segment, AL, of the coolant channel by 

Z 4AL/2 

c 



^B2 - %1 " "^-^ j "is^Z) dZ ^^^^ 



where Z is the center of the increment under consideration. Defining an 
c 

average surface heat flux for the increment, q^' = ^s^^s' incremental 
bulk temperature rise becomes 

ttD AL „ 

^B2"^Bl=S-r~%- (^8) 
P 

The local surface heat flux, <lg(Z), i ' related to the local vc>lumetric 
thermal source strength by observing chy'" all of the energy generated in the 
fuel is transferred out through he suriace of the cladding; i.e.. 
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I'D, dZ q^'(Z) - J d| dZ q"'(2) 



(49) 

This 3implifies to 



1 1 Df III 



(Z) q"'(2). (50) 



Similarly , 

^ _iii 



s 4D" • (51) 

s 
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4.8 Assessment of Local Coolant Conditions 

The current generation LMFBR's are designed to operate with no boiling 
of the liquid metal coolant. This requires that at every point In the flow 
the liquid metal temperature must be lower than Its saturation temperature 
at the local pressure. 

Tne maximum coolant temperature at any cross section of the reactor 
flow channel occurs In the fluid layer Immediately adjacent to the cladding. 
The magnitude of this temperature approaches the outer surface temperature 
of the cladding. Thus, the criterion for the absence of boiling in the reactor 
Is that the cladding surface temperature at all points be lower than the 
saturation temperature of the liquid metal at that level In the reactor. 

The axial location of the maximum cladding surface temperature can be 
obtained by writing an expression for T^CZ) in terms of q'^'* D^, D^, 
m, Cp, and h, differentiating with respect to Z, setting the derivative 
equal to zero, and solving for Z^. The maximum cladding surface temperature 
can then be obtained by evaluating T^Cz) at Z^. The details of this 
calculation are left as an exercise for the student. Unless the rate of 
decrease of the saturation temperature with axial position (due to frictional 
pressure drop) is appreciable, the critical point to check for surface boiling 

is near Z . 

m 

If the variation of cladding surface temperature and saturation 
temperature with axial position are known, it is a simple matter to visually 
check these data for the possibility of surface boiling. The output of the 
LMFBR Thermal- Hydraulics Code includes this information. 
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5.0 The LMFBR Thermal-Hydraulics Code 

The LMFBR Thermal-Hydraulics Code calculates the thermal-hydraulic 
performance parameters discussed in Section 4 for the LMFBR model described 
in Section 4. The basic operational steps of the code are listed below. 

1. Real<-* integer conversions 

2. Definition of statement functions 

3. Accept and print input 

4. Calculate flow area and equivalent diameter of 
coolant channel 

5. Determine calculation increment 

6. Calculate inlet pressure loss 

7. Initialize to 1st calculation increment 

8. Calculate bulk temperature rise and average temperature 
for 1st increment 

9. Calculate average coolant properties for 1st increment 

10. Calculate pressure loss for 1st increment 

11. Calculate average pressure and saturation temperature 
for 1st increment 

12. Calculate heat transfer coefficient for 1st increment 

13. Calculate cladding outer surface temperature for 1st 
increment 

14 • Check for surface boiling for 1st increment 

15. Calculate temperature distribution in fuel for 1st 
increment 

16 • Print 1st increment output 

17. Repeat steps 7-16 for remaining increments 

18. Calculate axit pressure loss and total pressure loss 
19 • Print coolant exit conditions. 
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The details of these specific steps requiring further explanation are 
found in the discussion below. A flow chart showing the calculations and 
logic is given in Section 5.4. 

Most of the statement functions are used for convenience in calculating 
property data. These statement functions are in the form of some property, 
either of the coolant, cladding, or fuel, as a function of the correspond- 
ing temperature. All of the function forms are polynomials obtained by 
fitting polynomials of various degree to tabulated data. The order of these 
polynomials varies, having been selected to insure that the fit is accurate 
to within 2% over the range of interest of the parameter. 

A detailed discussion of the input requirements is deferred until 
Section 5.2. For the moment it suffices to note that this input is 
supplied on three cards; one each for the geometry, the inlet flow conditions, 
and the reactor power level. 

Since the cross section of the model flow channel does not change with 
axial position, the flow area and equivalent diameter are constant axially. 
The flow area is calculated from 

A^ = f S2 cos 30--fD^-iD^. (3) 



The wire wrap diameter is related to S and D by 



D = S - D . 
WW s 



(52) 
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The equivalent diameter is then calculated from 

4A 



D « — - (5) 
e p 

where the wetted perimeter is given by 

In the code, the model coolant channel and fuel rod are sliced into 
short segments which are stacked axially to form the proper length core. 
The actual length of these segments is selected as a compromise between 
the very short segments for which the assumption of constant properties over 
the section axially is accurate and the very long segment which minimizes 
the computer time. An increment size equal to 1/100 of the active core 
length is built into the code. To study the effect of increment size, or to 
save computer time, this can readily be changed by altering two lines of 
the code. 

After the increment size has been computed, the coolant conditions at 
the beginning of the first increment are established. The bulk temperature 
and mass velocity at the entrance of the first increment are set equal to 
the core inlet conditions provided in ^he input. The pressure at the en- 
trance of the first calculation increment is obtained by subtracting the 
inlet pressure loss, computed from Equation 29 and Equation 37, from 
the pressure at the inlet of the core. In addition, counters required in 
the code logic are set equal to 1, denoting the first increment. 

The first major calculation in the code is the determination of the 
bulk temperature rise for the first calculation increment. This is calculated 
from Equation 48 with q^ determined from Equation 51. For small calculation 
increments, the averagjB volumetric thermal source strength can be accurately 
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approximated by evaluating Equation 1 at -H/2 + AL/2; i.e., at the center 
of the increment. 

The average bulk temperature, T^, for the increment is taken as the 
arithmetic average of the inlet and exit values. All coolant properties 
used in the pressure drop and heat transfer coefficient calculations are 
assumed to be functions of temperature only and are evaluated at the average 
bulk temperature. 

Both the friction factor and heat transfer coefficient calculations 

require the Reynolds number of the flow which is calculated from Equation 

33. Then f is calculated from Equation 34. To this value the 
circ 

multiplicative correction factors given by Equations 35 and 36 are applied 
to yield f for the actual reactor flow. Knowing f, the pressure drop for 
the increment is calculated from Equation 29. The Nusselt number is cal- 
culated from Equation 42 and the heat transfer coefficient is then backed 
out of the Nu definition. Equation 39. 

The average outer surface temperature of the cladding for the increment 
is determined from Newton's law of cooling with q^, T^ and h available from 
previous calculations for the increment; i.e.. 



II 



X = f + ^ . (53) 
s D n 

At this point, the code checks whether boiling is present in the increment. 
If boiling is observed an informative message is printed and the calculations 
are terminated. 

After the coolant channel calculations for the increment have been 
performed, the radial temperature distribution in the fuel is calculated for 
the increment. 
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First the temperature at the outer surface of the fuel is found. In 
the fuel rod model'^ the outer surface of the fuel is at the same temperature 
as the inner surface of the cladding; therefore the outer surface fuel 
temperature, T^, can be obtained by adding the temperature drop across the. 
cladding to the outer surface temperature of the cladding, f^, which was 
calculated earlier. 

The temperature drop across the cladding is obtained from Equation 25. 
In Equation 25 the thermal conductivity of the cladding appears- This 
property is temperature dependent and should be evaluated at the mean tem- 
perature of the cladding- However the mean cladding temperature is dependent 
on the temperature rise across the cladding- Thus, an iterative-type of 
solution is called for- The procedure used is to evaluate k^_^ at T^ , 

calculate T^, calculate f - , , evaluate at T - , compare 

<f f . clad mean c-II clad mean 

k to k^^i* ^-Plf k ^^>is within 1% of k ^, accept value of T^- If 
c-II c-I ^""^l ^ ^ 

k is not within 1% of k ^, let k ^ " ^ tt repeat the sequence - 
c-II c-I c-I c-II 

In this way, a cladding thermal conductivity accurate to within 1% is 
obtained- 

The next step is to calculate the temperature at the center of the fuel 
from Equation 17- Again a temperature dependent thermal conductivity, k^., 
appears, and an iterative calculation procedure like that described in the 
previous paragraph is' used- The fuel temperatures at the quarter, half, 
and three-quarter radii are then calculated from Equation 15- 

After the calculations for the first increment are completed the results 
are printed out- Discussion of the output is deferred until Section 5-3. 
It should be noted here, however, that in the output tabulation the tempera- 
tures are reported as having occurred at the axial center of the increment- / 
In other words, for the first increment the results are reported at 
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Z = -H/2 + AL/2. After the results for the first increment have been 
printed out the code increments to the next segment, and continues this 
process until all 100 segments have been spanned. 

After the completion of the calculations and printout for the 100 
segments of the core, the core exit conditions are calculated. This 
requires a calculation of the exit pressure loss using Equations 29 and 
37. Then P and T are computed and printed, eqding the computations. 
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5.1 Nomenclature 



Analysis Sytnbol 



Code Symbol 



Description 



^1' S' ^4 



P 
D 



WW 



e 



o ci. 



rc 



circ 



AFL 



CP 

DEQ 
DFU 
DCL 
DWW 
DR 



ERR 
F 
FRAT 

FCIRC 
FWWM 



Area 

Flow area of coolant channel 

Heat transfer area normal to 
n-direction 

Heat transfer area at radius r 

Heat transfer area at outside of 
cladding 

Integration constants 

Constant pressure specific heat 

Diameter 

Equivalent diameter of coolant channel 

Outside diameter of nuclear fuel 

Outside diameter of cladding 

Wire wrap diameter 

Ratio of wire wrap diameter to 
equivalent diameter 

Energy release per fission 

Error ratio 

Friction factor 

Ratio of unwrapped fuel bundle to 
circular tube friction factor 

Friction factor for circular tubes 

Ratio of actual to unwrapped fuel 
bundle friction factor 



Dimensional constant =* 32.174 



Ibm-ft 
2 

sec -Ibf 



EKLC 



Analysis Symbol 
G 

h 
H 

HP 
1 

k 

k 

c 
c-I 

^ TT 



K 

L 

e 
m 

N 



ff 

Nu 
Option 
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Code Symbol 

GIN 

HTC 
H 
HE 
HP 

ENTH 



KBC 
KCL 
KCLl 

KCL2 

KFU 
KFUl 

KFU2 



DELL 
MD0T 



NU 
0P 



Description 
Mass velocity 

Convectlve heat transfer coefficient 
Height of active core 
Extrapolated height of core 
Helical wire wrap pitch 
Enthalpy 

Enthalpy at Inlet of core 

Thermal conductivity 

Thermal conductivity of the coolant 

Thermal conductivity of cladding 

Thermal conductivity of cladding - 
at pre-lteratlon average cladding 
temperature 

Thermal conductivity of cladding - 
at updated average cladding 
temperature 

Thermal conductivity of fuel 

Thermal conductivity of fuel - at 
pre-lteratlon average fuel temperatur 

Thermal conductivity of fuel - at 
updated average fuel temperature 

Flow resistance coefficient 

Extrapolation length 

Length Increment 

Mass flow rate 

Number of fissionable fuel nuclei 
per unit volume 

Nusselt number 

Option « 1 for cosine, Option - 2 
^^for coupled program 



I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 

J 
I 



I 



I 
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Analysis Symbol Code Symbol Description 

p PW Wetted perimeter 

P PAV Pressure 

Pi Pressure at beginning of Increment 

P^ P2 Pressure at end of Increment 

P PEX Core exit pressure 
ex 

P, PIN Core Inlet pressure 
In 

Pe PE Peclet number 

Pr PR Prandtl number 

AP Pressure drop 

AP„ DELPE Exit or entrance pressure loss 
£ 

AP„ DELPF Frlctlonal pressure loss 
F 

q Heat transfer rate 

Heat transfer rate at 

q^ Heat transfer rate In n-dlrectlon 

a Heat transfer rate at D 

q" SHFAV Heat flux at 

s 

q''' QTPAV Volumetric thermal source strength 

q" QTPO Volumetric thermal source strength 

° at center of core 

r Fuel pin radial coordinate 

r^ Radius of fuel 

r Outer radius of cladding 
s 

R Core radial coordinate 

Re RE Reynolds number 

S S Center-to-center spacing between fuel 

rods 
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Analysis Symbol 



Code Symbol 



Description 



S/D 



^1/4 
r 

^3/4 



Bl 



B2 



ex 



in 



sat 
T 

w 
U 

"b 

X 

n 



m 



SRAT 

T 

TO 
T14 
T12 
T34 
TBAV 
TBI 

TB2 
TF 
TEX 
TIN 
TSAV 
TSATAV 



UB 



2C 



Ratio of fuel rod spacing to fuel 
rod diameter 

Temperature 

Temperature at center of fuel rod 

Temperature in fuel at r « 1/4 r^ 

Temperature in fuel at r = 1/2 r^ 

Temperature in fuel at r « 3/4 r^ 

Coolant bulk temperature 

Bulk temperature at beginning of 
increment 

Bulk temperature at end of increment 

Temperature in fuel at 

Core exit temperature 

Core inlet temperature 

Cladding surface temperature 

Saturation temperature 

Wall temperature 

Velocity 

Average velocity in flow channel 
Coordinate in n-direction 
Axial coordinate 

Axial location of center of increment 

Axial location of maximum cladding 
surface temperature 

Roughness size 

Angular coordinate 
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Analysis Symbol Code SyiTibol Description 

p VXS Absolute viscosity 

V Kinematic viscosity 

p RH0 Density 

5 Average fission microscopic cross 

f section 

T Time 

Neutron flux 

superscript - Average for calculation increment 
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5.2 Code Input 

The input to the LMFBR Thermal-Hydraulics Code is intentionally 
extremely simple. Only three cards are required; the first lists the 
geometrical quantities, the second lists the coolant inlet conditions, 
and the third lists the parameters that describe the volumetric thermal 
source strength distribution in the reactor. 

The required input data cards are listed below. The units of each 
parameter and their input format are included. Sample input data cards 
are shown in Figure 2. 

Card 1 — Geometry 

(in.) F 10.4 

D (in.) F 10.4 

s 

S (in.) F 10.4 

H? (in.) F 10.2 

H (in.) F 10.1 

Card 2 — Coolant inlet conditions 

G (Ibm/hr-ft^) E 10.3 

T^^ (F) F 10.1 

P. • (psia) F 10.1 

m 

Card 3 — Power distribution 

q' " (Btu/hr-ft*^) E 10.3 
o 

H (in.) F 10.1 

e 

Option (pure niamber) I 10 

On Card 3 the axial volumetric thermal source strength distribution 
through the core is specified. For independent operation of the code 
(Option «= 1) a cosine distribution is built into the code and q'' ' and H 
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Card 3 



Card 2 



Card 1 



0 ii[> n nnmm ooomioooQQjiiOQiiooiKioiiysiiooooDoyoDDooiiooooociioycDODOiicodH 

u H i n Mi0iiiniiri:iiiiiii;tiiiiiK<a?iniiaM jii !iiMri«iiiiiiOtiOowoiii/a<>»iiuiisM«i;M!»itiH UHii»ii 



r 



).40oe Oi? 



165.0 



IJI!0j00JQ0ll00l!HOj00ODC0Htl/i^!l!OnOOO 0 0 0 OC00OQIIO0CO00IlllO 0 0 0 ll3y0y0CO 0 0 0 0 llODD3 
lilts nH:giiiiniM!iiiiitiKi>i)nninNn!iN]Dii:H))4is!(]!]iii(c«ii!u«ouiiu<)!irjis;s)s*!i»iiM!iHiii;u;uitiiTUt)^ 



O.SSOd 0,£5i)0 0,3250 



15.00 

0 



60. i 



0 0 0 ^000^,'J(!OJOIlHHIlQjSOOOJ]|)OOOOIljiOOQOQI10j;i)DII()OIiDIIDri(!HGI!llOOOOOI)llO[!90bDO 
n 1 u I ) I iioiiin)in)iii)iiiiiir»:i2i».H!i»N]oii.ij]ios»njt]i4ii<ioou*M;t)u(jiosis!uinist!iMiut(ii)ii(uiitnti» 



' 1 1; 



2 222222cc22222222c2222222222c222^2222c2222222222222n222222?3222222n2222?2222222 , 

0 u 

0 3 3 3 5 3 , 3 3 3 3 3 3 3 3 3 , 3 3 3 3 3 3 3 3 3 ,, 3 333?33333i333 3 3 3 3 3 3 37 3 3 3 3 33n3 3 3 3 3 nJ 3 3 3 3 ijnn 3 3 3 3 } 



3 44UUU444444U4H44U444U4MUU44U4^4M;n44MU4U4U44UUUMUUUU44< ^ 
s c 

8 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 a!i 5 5 5 5 !i^5 5 5 o5ll 5 5 5 5 5 5 ^5555S 5 5 5 5 5 i 5 5 5 5 5 5 5 5 5 5 5 5 5 S55!i55S 5 5 5 !ij:i5^ ^ 

\\\\\ wwmw seuiS888 (ussissisc sitssinis nssnsssutssetssississssssjsss 



g5 Figure 2. LMFBR Ttiemial Hydraulics Code Input. 
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are the parameters which indicate its level and period- When the LMFBR 
Thermal-Hydraulics Code is coupled with other codes (Option = 2) the 
cosine distribution is over-ridden by the actual power distribution which 
is supplied by the main program. 
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5.3 Code Output 

Three sets of information are printed out. Each is printed under 
completely explicit headings sucli that there should be no interpretation 
difficulties. 

The first set of output information consists of a listing of the input 
that was supplied to the code. 

The second and main set of output is a tabulation of the temperatures, 
h, and AP for each increment in the core. The temperature printed out are 

'^o* '^1/4* '^1/2' ^3/4' '^f '^s* '^B' '^sat* ^^o"^ this data, observations of the 
maximum fuel and cladding temperatures and their locations can be made. 

The third set of data reports the core exit pressure and temperature. 
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5.4 Code Flow Chart 



& 



kf(T), k^(T) 
kgCT), p(T) 
P(T), Cp(T) 



Input: 
Df, Dg, S, HP, H 

^q^". Hg, Option 



OuLput: 
Df, Dg, S. HP, H 

^ » ^In' ^in 
, q^". H ^. Option 



cos 30 



WW' 



p = i^D + D ) 
*^w s ww^ 



Functional relations between 
variables supplied to code 
as polynomials. 



Geometry, inlet coolant 
conditions, and power 
distribution supplied 
to code. 



Printout of geometry, 
inlet coolant conditions, 
and power distribution. 



Calculation of flow channel 
geometrical parameters. 



100 



Selection of calculation increment. 
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U3=G /p(T^^) 



Re = 



f , = 0.184/Re 
clrc 



0.2 



o clrc f 



ik) 



circ s 



o e 



AP = 23f — ilL-B. 

F "^-^^ D 2g 
e c 



P = - AP^ - AP^ 
1 In E F 



= 0.5 



Calculation of permanent pressure 
loss at inlet and resulting pressure 
at start of active core. 



^Bl ^in 



Z = -H/2 + AL/2 
c 

1 = 1 






nz 






q' " = q" ' 


cos^j-S- 




e 


T = T + 
B2 ^Bl 


i'" nD^/8mCp(Tgj^) 


T = (T + 
B Bl 


V/2 




Initializing temperature, 
position, and increment 
counter for 1st increment. 



Calculation of average bulk 
temperature of coolant for 
increment. 
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f ^ = 0.184/Re 
clrc 



0.2 



f = f ^ T 
o clrc f 



circ ^ s 



o e 



25 



AP = f AL P<Jb2"b 

^ D 2g 
e '^c 



P = P - At 

T = T (P) 
sat sat 



Pr = 



^B.(^B> 



Pe = RePr 



Nu = Nu(Pe, — ) 

- Nuk^T^, 
h » 



T = T + 

s B 4 d;; h 

8- 



Calculation of pressure drop 

and average saturation temperature 

for Increment. 



Calculation of heat transfer 
coefficient and average cladding 
surface temperature for Increment. 
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1 — ^ 
— ^at 




>■ 


k = k (T ) 
cl c s 








Output; 
''Surface 
Boiling" 





q"'Dj D 

^s-^ 8k ; 

cl f 


f = 

c 


(Tg + T^)/2 










^1 - ^c2 




e = 


^c2 





Check to see if subcooled boiling 
is possible. If possibility is 
detected, a message is printed and 
calculations are terminated. If 
no boiling iR indicated, program 
proceeds to calculation of fuel 
and cladding temperatures. 



Calculation of temperature rise 
across cladding. The thermal 
conductivity of the cladding is 
evaluated at the mean cladding 
temperature by an iterative 
techn j que . 




^fi = ^f(^f> 




T = T + — 

o f 16k^j^ 



T. = (T. + T )/2 
t to 



k,2 = (T,) 




k = k - 
cl c2 



Calculation of fuel centerllne 
temperature. The thermal con- 
ductivity of the fuel is eval- 
uated at a mean fuel temperature 
by an iterative technique. 



ERIC. 
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k - k 


e = 












^1/4 


« T 

o 


256 








^1/2 


= T 

o 


64 k,^ 








^3/4 


= f 

o 


256 k^^ 




Evaluation of fuel temperature 
at the quarter, half, and three 
quarter diameter points. 



Output: 



^I, Z^, T^, T^/^, ^3/4* 

^s^ ^B^ ^Sat» 

AP, CHFR / 




1 = 1 + 1 



Z = Z + AL 
c c 

^Bl = ^B2 

^ = ^2 



0 



■0 



Printout of thermal-hydraulic 
parameters for increment. 



Check to see if entire core has 
been calculated. 

Initializing for next calculation 
increment. 
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T = T_, 
ex B2 



B p(T ) 
ex 



Re = 



U„ D p(T ) 
B e ^ ex 



f = 0.184/Re 
clrc 



0.2 



o clrc f 



o /S 



clrc \ s 



o e 



AT P(T ) 



P<^ex) "b 
2g. 



P = P. - AP_ - AP^ 
ex 2 F E 




Calculation of pressure loss In 
gas plenum and downstream rod 
support region. 



Printout of coolant core exit 
conditions . 



ERIC 
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5.5 Code Examples 

The following examples demonstrate the use of the LMFBR Thermal- 
Hydraulics Code to calculate the thermal-hydraulic behavior of a 
typical LMFBR. 
Example 7. 

A LMFBR core is composed of 0.250 in. diameter fuel rods 5 ft long 
and spaced 0.325 in. in a triangular array. The fuel diameter is 0.220 in. 
and the wire wrap pitch is 12 in. The liquid sodium inlet conditions are 
800 F and 200 psia. The average coolant velocity in the Inlet of the core 
is 28.80 ft/sec. The average fuel rod produces 8.04 kw/ft. Assume an 
extrapolation length equal to 5% of the core height. Determine the mag- 
nitudes and locations of the maximum fuel temperature and the maximum 
cladding temperature. 
Solution. 

Xhe mass velocity, extrapolated core height, and peak volumetric 
thermal source strength must first be computed from the data given, 

G = PU^ = (53.1 Ibm/ft^) (28.80 ft/sec) (3600 sec/hr) 
B 

= 5.5 X 10^ Ibm/hr-ft^ 

H =H+2L=H+ 2(0. 05H) = I.IOH 
e e 

= 1.10 (5 ft) = 5.5 ft = 66.0 in. 

q'"=2P /[H sln(nH/2 H ) ] 
^ o ave e f e 
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2(8.04 kw/ft)(12 ft) (3413 Btu/kw-hr) /{ (5.5 ft) (0.220 in.)^ 



(ft/12 In.)^ sin[]l(5 ft)/2(5.5 ft)]} 



1.5 X 10^ Btu/hr-ft'^. 



From the output printout of the code (see next five pages) we find 

T , , = 1120.0 F at Z = 1.925 ft 
clad max 

T. . = 3548.2 F at Z = 0.025 ft , 
fuel max 
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LMFB^ THERMAL HYDRAULICS CODE 

* « 4: * « « « « « 4t « « « 4< * A 4c >: ■' « ♦ # « « « 4i 4> « 

INPUT UATj!.. 



FUEL 0.0. 
( INCH 1 
0.2? CO 



CLAr. n.D, 

( INCH ) 
0.2300 



ROD PITCH 
( INCH ) 
0.3250 



ACTIVE CORE LENGTH 
( INCH ) 
60.00 



WIRE WRA** PITCH 
( INCH ) 
12.00 

INLET M^SS VELOCITY 
( L3/HR-FT**2 ) 
O.5500E 07 

VOLUMETRIC THERMAL 
SOURCE STRENGTH 
( BTU/HR-FT**3 ) 
0.15005 09 



CORE INLET TEMPERATURE 
( F I 
800.00 

EXTRAPOLATED HEIGHT 

(INCH ) 
66.00 



CORE INLET PRESSURE 
( PSIA ) 
200.00 



OPT I JN 



EN.) OF INPUT DATA 
41 « « 4[ « 4r 4i 1c 4i 4c « « « 4e « # « * « >f « « * >C * * * * * 



AT THE INLET CF THE CORE 



CHANNEL AVE. VELOCITY REYNOLDS FRICTION FACTOR 

( FT/SEC) 

28.80 111627.4 0.0180 
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NO Z-LOC. 10 Tl/4 Tl/2 T3/4 TF 

T CLAD T BULK T SAT H T CDFF PCLTA P 

1 -2.%75 1130. i 1058.1 966.1 813.3 

tiCZ.b 300.3 2289.8 44078.8 50.31 

2 -2.425 1^12.1 1201.0 1113.2 989.5 816.3 

y03.6 801.1 2289.2 44068.6 50.31 

3 -2.J75 I2rip,.3 1275.1 1171.1 1024.6 819.5 

804. b (jJl.9 2288.6 44056.6 50.30 

4 -2.J25 136B.3 1353.0 1231.9 1061.4 822.7 

^^06.^? 802.9 2268.0 44042.9 50.30 

5 -2.275 1452.1 1434.5 1295.6 1100.0 826.1 

tiJ7.6 803.9 2287.4 44027.6 50.30 

6 -2.225 1341.2 1521.2 1363.3 1140.9 829.5 

809.2 805.1 2286.8 4401C.7 50.30 

7 -2.175 lo33.1 1610.7 1433.1 1183.1 833.1 

HIO.'.. 800.4 2286.2 43992.1 50.29 

8 -2.125 1728.4 1703.4 1505.5 1226.8 836.7 

dl2.7 807.8 2285.6 43971.9 50.29 

9 -2.075 It526.5 1798.8 1579.9 1271.8 840.4 

:JU.6 809.4 2285.0 43950.0 50.28 

10 -2.025 l-J2b.i 1895.9 1655.7 1317.6 844.2 

nio.t eil.O 2284.4 43926.6 50.28 

11 -1.975 20<'6.7 1993<.6 1732.1 1363.7 848.0 

SlH.l 812.8 2283.8 43901.7 50.28 

12 -1.925 2126,5 2090.7 1807.9 1409.6 852.0 

820. bl4.7 2283.2 43875.2 50.27 

13 -1.J7S 2224.3 2185.9 1882.2 1454.7 856.0 

823.3 bl6«6 2282.6 43847.2 50.27 

14 -1.625 2318.9 2273.0 1954.2 1498.3 860.1 

325.7 618.7 2282.0 43817.6 50.26 

15 -1.775 2'»06.6 2363.3 2021.0 1539.0 864.3 

d26.2 820. V 2281.3 43786.7 50.26 

16 -1.725 2493.4 2447.8 2087.2 1579.4 868.5 

330.9 823.2 2280.7 43754.3 50.25 

17 -1.675 2575. C 2527.2 2149.4 1617.5 872.8 

833.6 825.6 2280.1 43720.4 50.25 

18 -1.625 2651.0 2601.2 2207.5 1653.2 877.1 

336.4 628.1 2279.5 43685.2 50.24 

19 -1.375 2722.1 2670.4 2261.9 1686.8 881.?? 

839.3 830.7 2278.9 43b4e.7 50.23 

20 -1.525 2787.3 2733.9 2312. C 1717.8 886.0 

842.? 333.4 2278.3 43ol0.8 50.23 

21 -1.475 2d<^8.5 2793.5 2359. u 1747.1 890.5 

845.4 836.2 2277.7 43571.7 50,22 

22 -1.425 2905.3 2849.3 2403.1 1774.7 895.0 

y4i3.C 839.1 2277.0 43531.3 50«21 

23 -1.375 29t.2.4 2904.5 2446.7 1P02.1 899.6 

{,'51.9 842.0 2276.4 43489.7 50.21 
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2^ -I. 32b 3011.2 2952.0 2484«4 1826.0 904.2 

8bb.4L 845.1 2275.8 43446.9 50.20 

25 -I.27*J -.iObb.l 2995.8 2519.3 1848.3 908.9 

Sbti.f- 848.2 2275.2 43403.0 50.19 

26 -1.225 3098.5 3037.1 2552.2 1869.4 913.5 

862.1 B51.5 2274.6 43357.9 50.19 

27 -1.175 31t2.9 3080,4 2586.7 1891.5 918.2 

86ti.7 854.8 2273.9 43311.9 50.18 

28 -1.125 iie2.2 3118.7 2617.4 1911.4 923.0 

869.3 858.1 2273.3 43264.8 50.17 

29 -1.075 3213.2 3149.0 2641.8 1927.6 927.7 

H73.0 861.6 2272.7 43216.7 50.17 

30 -1.J25 J2«i9.5 3184.5 2670.3 1946.2 932.5 

b76.b 865.1 2272.1 43167.7 50.16 

31 -0.975 32B1.3 3215.4 2695.3 1962.8 937.3 

B8J.6 868.7 2271.5 43117.8 50.15 

32 -0.925 3306.3 3239.9 2715.3 1976.4 942.1 

884. i> 672.4 2270.8 43067.0 50.15 

33 -0.b75 33^2.4 3265.4 2736.0 1990.5 946.8 

809. b 876.1 2270.2 43015.4 50.14 

34 -0.825 3361.6 . 3293.9 2759.1 2006.0 951,6 

892.4 879.9 2269.6 42963.1 50.13 

35 -0.775 33B1.J 3313.2 2775.1 2017.3 956,4 

896.5 883.7 2268.9 42910.0 50.12 

36 -0.725 3-^02.0 3333.5 2791.8 2029.1 961,2 

900.6 887.6 2268.3 42856.3 50.12 

37 -0.675 J425.6 3356.5 2810.7 2042.1 966.0 

904.7 891.6 2267.7 42802.0 50.11 

38 -0.625 3444.2 3374.8 2825,9 2052.9 970.7 

9U8.Q 895.6 2267.1 42747.1 50.10 

39 -0.575 3456.6 3386.9 2836.3 2061.0 975.5 

913.1 899.6 2266.4 42691.6 50.10 

40 -0.52b 3471.1 3401.1 2848.4 2070.0 980.2 

917.3 903.7 2265.8 42635.7 50.09 

41 -0.475 3489.1 3418.7 2863.0 2080.5 984.9 

921.6 907.9 2265.1 42579.3 50,08 

42 -3.425 35C1.7 3431.1 2873.7 2088.6 989.6 

925. b 912.0 2264.5 42522.5 50.08 

43 -0.37b 3PC8.1 3437.5 2879.6 2094.0 994,2 

933.1 916.2 2263.9 42465.4 50.07 

44 -0.325 3517.2 3446.5 2887.6 2100.6 998.8 

934.5 920.5 2263.2 42408.0 50.07 

45 -3.275 3529.5 3458.6 2898.0 2108.6 1003.4 

938. y 924.7 2262.6 42350.3 50.06 

46 -J. 22b 3536.6 3465.6 2904.4 2114.2 1007.9 

943.2 929,0 2262o0 42292.5 50.06 

47 -0.175 3j42.4 3471.3 2909.9 2119.3 1012.4 

947.5 933.3 2261.3 42234.4 50.05 

48 -0.12b 3=342.7 3471.8 2911.2 2121.9 1016.8 

951. 937.6 2260.7 42176.3 50.04 

49 -0.37b 3545.7 3474.8 2914.6 2125.7 1021.2 

956.2 941.9 2260.0 42118.2 50.04 

50 -0.025 35A7.6 3476.7 2917.0 2128.9 1025.5 

960.6 946.2 2259.4 42060.0 50.04 
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51 0.025 3t>48.2 3A77.5 2918.6 2131.6 1029.8 

9(>4.V 950.5 2253.8 42001.9 50.03 

52 0.075 3547.7 3477.1 2919.3 2133.7 1034.0 

9iW.-:) 954,9 225d.l 41943.9 50.03 

53 0.125 3546.1 3475.6 2919.1 2135.3 1038.1 

973. t 959.2 2257.5 41386.0 50.02 

54 0.175 3543.2 3473.0 2918.0 2136.4 1042.2 

977.9 963.5 2256.8 41828.3 50.02 

55 0.225 3543.1 3473.0 2918.9 2138.6 1046.2 

932. 967 8 2256.2 ^1770. 8 50.02 

56 0.275 3537.6 3467.7 2915.7 2138.4 1050.1 

9B6.5 972.1 2255.5 41713.6 50.01 

57 0.325 3530.9 3461.3 2911.7 2137.6 1053.9 

990.7 976.4 2254.9 41656.7 50.01 

58 0. 375 3ijl9.0 3449.9 2903.7 2134.5 1057.7 

994.9 9d'J.7 2254.2 41600.2 50.01 

59 0.425 3510.4 3441.6 2898.1 2132.8 7.061.4 

999.1 934.9 2253.6 41544.1 50.01 

60 0.4^5 35J0.5 3432.1 2891.6 2130.5 106/:.9 

1003.,? 989.1 2252.9 41488.5 50.00 

61 0.525 3492.9 3424.8 2886.8 2129.2 1068.4 

10U7.3 993.3 2252.2 41433.4 50.00 

62 0.575 3479.9 3412.3 2877.9 2125.4 1071.8 

1011.3 997.5 2251.6 41378.3 50.00 

63 0.625 3462.0 3395.0 2865.3 2119.4 1075.1 

1015.3 1001.6 2250.9 41324.8 50.00 

64 0.675 3447.2 3380.6 2855.0 2114.7 1078.4 

101';.? 1005.6 2250.3 41271.5 50.00 

65 0.725 3434.6 3368.5 2846.3 2170.9 1081.5 

1023.1 1009.7 2249.6 41218.8 50.00 

66 0.775 3416.2 3350.7 2833.2 2104.6 1084.5 

102O.9 1013.6 2249.0 41166.8 50.00 

67 0.325 3392. R 3328.1 2816.4 2096.0 1087.3 

1030.7 1017.6 2248.3 41115.6 50.00 

68 0.375 3372.5 3308.4 2801.9 2088.7 1090.1 

1034.3 1021.4 2247.6 41065.2 50.00 

69 0.925 3^54.2 3290.7 2788.9 2082.2 1092.8 

103t!.L 1025.2 2247.0 41015.6 50.00 

70 0.975 3329.9 3267.1 2771.2 2073.0 1095.4 

1041.5 1029.0 2246.3 40966.9 50.00 

71 1.025 3300.4 3238.5 2749.7 2061.4 1097.8 

1045.0 1032.7 2245.6 40919.2 50.00 

72 1.075 3273.9 3212.8 2730.5 2051.2 1100.1 

1048.4 1036.3 2245.0 40872.3 50.00 

73 1.125 3245.6 3185.4 2709.8 2040.0 1102.3 

1051.7 1039.8 2244.3 40826.4 50.00 

74 1.175 3218.1 3158.7 2689.7 2029.2 1104.4 

1054.9 1043.3 2243.6 40781.6 50.00 

75 1.225 31b5.1 3126.7 2665.4 2015.8 1106.4 

1058.1 1046.7 2243.0 40737.8 50.01 

76 1.275 3147.4 3090.2 2637.6 2000.4 1108.2 

1061.1 1050.0 2242.3 40695.1 50.01 

77 1.325 3111.6 3055.3 261J.2 1985.6 1109.9 

1064.1 1053.2 2241.6 40653.5 50.01 
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78 1.375 3076. 'j 3021.3 2585.3 1971.2 1111.5 

1066.9 1056.4 2240.9 43613.1 50.01 ■ 

79 1.425 3034.8 2980.8 2554.3 1953.8 1113. 0| 

1069.7 1059.5 2240.3 40573.8 5C.0I 

80 1.475 ^'9^0. a 2933.1 2521.7 1935.3 1114.3 

I07Z.4 1062.4 2239.6 <;0535.7 50.02 I 

81 1.525 2y''f2.0 2890.7 2485.4 1914.6 1115.51 

IC74.9 1065.3 2238.9 40498.9 50.02 

82 1.575 2-JS4.J 2844.1 2449.7 1894.2 1116. 

1077.4 1068.1 2238.:: 40463.3 50.02 

83 1.625 2846.0 2797,4 2413.9 1873.7 1117.4 

1079.8 1070.8 2237r5 40429.0 50.02 . 

84 1.675 2790.1 .743.2 2372.? 1849.7 IIIS.ZI 

10.02.0 1073.4 2236.9 40396.0 50.03 ■ 

85 1.725 2731.5 2686.2 2328.3 1824.4 1118.9 

1084.1 1075.9 ?.:^6.2 40364.4 50.03 I 

86 1.775 2663.0 2624.5 2280.8 1796.9 J119.4I 

1086.2 1078.2 2235.5 40334.0 50.03 

87 1.825 2606.9 2565.2 2235.1 1770.4 1119.7| 

1088.1 1080.5 2234.8 40305.2 50.03 | 

88 1.875 2536.4 2496.6 2182.3 1739.7 1120.0 

1C39.9 1082.7 2234.1 40277.6 50.04 - 

89 1.925 2464.0 2426.2 2128.0 1708.0 1120. 0| 

1091.6 1084.7 2233.4 40251.5 50.04 ■ 

90 1.975 23fa7n3 2351.7 2070«5 1674.5 1120.0 

1093.1 1086.7 2232.7 40226.9 50.04 I 

91 2.025 2316.9 2273.5 2010.1 1639.1 1119. 8| 

1094. ft 1088.5 2232ol 40203.7 50.04 

92 2.J75 2^:23.1 2192.1 1947.2 1602.3 1119. 5| 

1095.9 1090.2 2231.4 40181.9 50.05 | 

93 2.125 2139.1 2110. 5> 1884.1 1565.3 1119.0 

1097.1 1091. 8 2230.7 40161.7 50.05 

94 2.175 2050.9 2024.7 1817. ^ 1526.4 111B.4| 

1098.2 ^0^3. 3 2230.0 401^3.0 50.05 ■ 

95 2.225 1961.2 1937.5 1750.3 1486.7 1117.7 

1099.1 1094.7 22?.9.3 40125.8 50.05 J 

96 2.275 1870.8 1849.6 1632.3 1446.6 1116.fl| 

llC'j.O 1095.9 2228.6 40110.0 50.05 

97 2.325 17b0.b 1761.9 1614.4 1406.6 lll5.a| 

1100.7 1097.1 2227.9 40095.9 50.06 I 

98 2.375 1691.3 1675.1 1-47.1 1366.9 lllV.6'" 

1101.3 1098.1 22L7.2 40083.3 50.06 

99 2.425 1603.8 1590.0 1481.2 1327.9 1113.31 

1101.7 1098.9 2226.5 40072.2 50.06 I 

100 2.475 15J 8.0 1506.6 14t6.5 ir.89.6 1111.9 

riU2.G 1099.7 2225.8 40062.7 50.06 | 

EXIT PKtSSURE EXIT TEMPEf.ATlIRE 



114.0 1100. 0 



I 



PDWtK PEK Ur,'IT LENGTH IKW/FT) 8.04 | 
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Example 8« 

The effect of oversized cladding on LMFBR thermal-hydraulic per- 
formance is to be studied. Consider the LMFBR described in Example 7 but 
with an outer cladding diameter of 0.275 in. (instead of 0.250 in.). 
Compare the following thermal-hydraulic parameters for the "oversized 
cladding" and "normal" flow channels: 

a. Location and magnitude of maximum fuel temperature 

b. Location and magnitude of maximum cladding temperature 

c. Exit sodium temperature. 
Solution . 

Each flow channel in the fuel assembly will experience the same 
pressure drop. The "oversized cladding" flow channel will have a smaller 
equivalent diameter and hence will pass less flow for the same pressure 
loss (as shown in Example 5) . To assure approximately the same pressure 
loss as in the "normal" channels the mass velocity of the ''oversized 
cladding" channel mxjist be adjusted. The equivalent diameters are 
computed first. 



II' 



'Normal" channel: 




— S cos 30 = -^(D + D ) 
2 OS WW 



1/2(0.325)^ cos 30^ - f[ (0.250)^ + (0.075)^] 



0.0191 in.^ 



w 



7(D + D ) = T(0.250 + 0.075) 

/ S WW ^ 



0.511 in. 



J) = 4 A-/p 
e f *^w 
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- 4(0.0150) /O. 5 11 
« 0.1497 in. 



"Oversized cladding" channel: 



A, = I S2 cos 30° - f (D^ + D^) 

= 1/2(0.325)^ cos 30° - |l (0.275)^ + (0.050)^1 



0.0150 In.^ 



p - ?(D + D ) - f (0.275 + 0.050) 
'^w 2 s WW t. 



= 0.511 In. 

= /Pw 

= 4(. 0150)70.511 

= 0.1173 in. 
From the results of Example 5, 



G = G„(D „ /D^„)^^^ 
oc N eoc eN 



= 5.50 X 10^ (0.1173/0.1497)^/^ 



6 2 
= 4.69 X 10 Ibm/hr-ft . 



The input to the LMFBR Thermal-Hydraulics Code is the same as that 
for Example 7 except for the new value of G computed above and the over- 
sized D . From the code (printout next five pages) the following results 



s 

are obtained 
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« 4> 4t « # « A 1(1 ;X « « « « » « « 4: « * 9^ « « « « « )|c « « « « 

LMF8R THERMAL HYORAULICS CODE 

4c 4c :;e It ^ :^ « jgc lie 1^ <c 4c * « 

INPUT DATA 



FUEL J.t). CLAU U.n. ROD PITCH ACTIVE CORF LENGTH 

I INCH ) ( INCH ) ( INCH ) ( INCH I 

0.22C0 0.2750 0.3250 60.00 



HIRE ^RAP PITCH 
( INCH ) 
12.00 



INLET MASS VELOCITY 
( L8/HR-FT*»2 ) 
0.4690E 07 

VOLUMETRIC THERMAL 
SOURCE STRENGTH 
( BTU/HR-FT**3 ) 
0.I500E 09 



CURE INLET TEMPERATURE 
( F ) 
SCO. 00 

EXTRAPOLATED HEIGHT 

(INCH ) 
66.00 



CORE INLET PRESSURE 
( PSIA ) 
200.00 



OPTION 1 



ENU .JF INI'UT DATA 



AT THE INLET HF THF- CORE 



CHANNEL AVE. VELOCITY REYNOLDS # FRICTION FACTOR 
(FT/SEC) 

24.56 75497.9 0.0195 

4c 4( 4> 4c 4c 4: >tc ^ 4c « « « 4: 9k u 4c t- ^ V 4c 4: 4c 4c li. 4: 4> 4> # 4c 4c 4c 4< * 4> 4> 4> 4> 4> 4c 4> 4c 4c 4c 4c 4> <t 4c 4> 4> 4< 4< 4: 4< 4i 4c 4> 4c 4c 4c 4c 4i 4c * 
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NO 2-LOC. TO Tl/4 Tl/2 T3/^ TF 

T CLAD T HULK T SAT H T COEF DELTA P ■ 

1 -2.47i llbO.B 1141.0 1068.1 965.4 821.6m 

802. b 800.5 2290.6 38036.8 51.59 ■ 

2 -2.425 1225.6 1214.4 1125.8 1001. 1 826.4 

804.? 801.6 2289.9 38024.1 51.59 

3 -2.376 1305.0 1291.7 1186.6 1038.6 831. 4H 

80i>-r 802.8 2289.3 38009.4 51.59 ■ 

4 -2.325 1-388.5 1373.0 1250.5 1078.0 836.5 

^807.7 804.2 2288.7 37992.7 51.58 ■ 

5 -2.275 1476. 1 1458.2 1317.5 1119.3 841. sB 

80S. 7 805.8 2288.1 37973.8 51.58 

6 -2.225 1569.1 1548.8 1388.6 1163.0 a47.2« 

811.9 807.6 2287.5 37953.0 51,57 

7 -2.175 1665.1- 1642.3 1462.0 1208.2 852.8"* 

814.2 809.5 2286.9 37930.1 51.57 

8 -2.125 1764.4 1739.0 1537.9 1254.8 858. 4B 

81o'.7 811.6 2286.2 37905.3 51.56 ■ 

9 -2.075 1366.2 1838.0 1615.7 1302.5 864.2 

819.3 813.9 2285-6 37878.5 51.56 ■ 

10 -2.025 1969.2 1938.3 1694.4 1350.9 870»lJ| 

822.2 816.3 2285.0 37849.8 51.55 

11 -1.975 2J 72 . 2 2 038 . 6 1 7 73 . 2 1 399 . 4 8 7 6.1- 

825.2 818.9 2284.4 37819.1 51.55 ■ 

12 - 1.925 2173. 7 2137.4 1850.8 1447.2 882.2* 

828.3 821.7 2283.8 37786.6 51.54 

13 -1. 67b 7.272.2 2233.3 1926.2 1493.8 888. 4l 

831.6 824.6 2283.1 37752.2 51.53 I 

14 -1.325 2J64.5 2323.2 1997.0 1537.7 894.7 

835.a 827.7 2282.5 37716.0 51.52 ^ 

15 -1 .775 2455.1 2411.4 2066.6 1580.9 ^^''^^l 

d33.7 831.0 2281.9 37678.0 51.52 

16 -1.725 25^0. 3 2494.4 2132.1 1621.9 907*5 

842.4 834.4 2281.3 37638.3 51.51 i 

17 - 1.675 2616.8 2569.0 2191.2 1659.0 

340.3 837.9 2280.6 37596.8 51.50 

18 -1.025 2694.8 2645.0 225lc3 1696.9 '■■■?.0«im 

nSJ.A 841.6 2280.0 37553,6 51.49 ^ 

19 -1.575 " 27f2.2 2710.7 2303.5 1730.2 

"854. o 845.5 2279.4 37506 .9 51.48 

20 > -1.525 2326. 3 2773.2 2353.3 1762.0 9''..2| 

%5d.'^>- - 849.4 2278.8 37462.5 51.47 " 

21 -1.475 2887.9' 2833.2 2401.2 1792. c 94i.I 

8b3.A 853.6 2278.1 37414.6 51.46 I 

22 ^-l.'tZ5 r29'^i3.4 2887.3 2444,5 1820.9 941. 9i 

868. C 857.8 2277. i» 37365.2 51.46 

23 -1.375 2997.5 2940.1 2486.8 1848.5 954.9 



872.7 A62.2 2276.9 37314.3 51.45 



I 
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51 0.025 3569.4 3501.4 2964.3 2208.0 1149-1 

1039.3 1023.8 2258.8 35509.5 51.28 

52 0.J75 3 5*5.6 3497.9 2963.0 2209.8 1155.3 

1045.8 1030.3 2258.1 35439-8 51.28 

53 0.125 3567.3 3500.2 2966.2 2214.2 1161.4 

1052.2 1036-7 2257.4 35370-2 51-28 

54 0.175 3565.5 3498.2 2966-0 2216.6 1167-4 

1058-.' 1043.2 2256.8 35301-0 51-28 

55 0.225 3561.9 3494.8 2964-7 2218-3 1173-2 

1065.1 1049.6 2256-1 35232.1 51.29 

56 0.275 3553.5 3486-8 2959.8 2217-8 1179-0 

1071.5 1056.1 2255-4 35163.6 51.29 

57 0.325 3551.0 3484.5 2959.4 2219.9 1184-6 

1077.9 1062.5 2254.8 35095.5 51.29 

58 0.375 3543.9 3477.8 2955.5 2219.9 1190.1 

1084.2 1068.8 2254.1 35027.9 51.30 

59 0.425 ;o32.1 3466.5 2948.0 2217.8 1195-5 

1090.4 1075.2 2253.4, 34960.8 51.31 

60 0.475 3526.4 3461.1 2945.0 2218.3 1200.8 

1096.7 1081.5 2252.8 34894.4 51.31 

61 0-525 3516.1 3451.2 2938.6 2216.6 1205.9 

1102.8 1U87.7 2252.1 34828.6 51.32 

62 0.575 3504.5 3440.1 2931.1 2214.3 1210-9 

1108.9 1093.9 2251.4 34763.6 51.33 

63 0.625 3491.6 3427.7 2922.7 2211.4 1215.7 

1114.9 1100.1 2250.7 34699.3 51.34 

64 0.675 3477.5 3414.1 2913.2 2207.9 1220.4 

1120.9 1106.2 2250.1 34635.8 51.35 

65 0.725 3458.9 3396.2 2900.4 2202.3 1224.9 

1120.7 1112.2 2249.4 34573.1 51.36 

66 0.775 3442.9 3380.7 2889.5 2197.7 1229.3 

1132.5 1118.1 2248.7 34511.4 51.37 

67 0.825 3428.3 3366.6 2879.6 2193.7 1233.5 

' 1138.2 1124.0 2248.0 34^50.6 51.38 

68 0-375 34vl6.0 3345.1 2863.9 2186.2 1237.5 

1143.8 1129.8 2247.3 34390.8 51.39 

69 0-925 3388.7 3328.4 2851.9 2180.8 1241.4 

1149.3 1135.5 2246.7 34332.0 51.41 

70 0.975 ?^7.0 3307.4 2836.5 2173.4 1245.1 

1154.7 1141.1 2246.0 34274.3 51.42 

71 1.325 3344.0 3285.1 2820.2 2165.4 1248.6 

1160.0 1146.6 2245.3 34217.8 51.43 

72 1.075 3316.2 3258.2 2800.2 2155.1 1252.0 

1165.2 1152.1 2244.6 34162.4 51.45 

73 1.125 3293.5 3236.2 2783.9 2146.9 1255.2 

1170.3 1157.4 2243.9 34108.2 51.46 

74 1.175 3266.1 3209.7 2764.1 2136.6 1258.2 

1175.2 1162.6 2243.2 34055.2 51.48 

75 1.225 32o7.0 3181.5 2743.0 2125.5 1261.0 

118J.0 1167.7 2242.5 34003.5 51.49 

76 1.275 3206. 1 3151.5 2720.5 2113.4 1263.6 

1184.7 117^.6 2241.8 33953-2 51.51 

77 1.325 3173.3 3119.7 2696.5 2100.4 1266.0 

1189.3 1177.5 2241.1 33904.2 51-52 

77 



- 73 - 



78 1.375 3138.5 3086.0 2671.0 2086.5 1268.3 

1193.7 1182.2 22^0.5 33856.5 51.5't 

79 1.425 3101.8 3050.3 2643.9 2071.6 . 1270.3 

1193.0 1186.8 2239.8 33810.4 51.55 

80 1.475 3062.9 3012.6 2615.2 2055.6 1272.2 

1202.1 1191.3 2239.1 33765.6 51.57 

81 1.525 3021.9 2972.8 2584.8 2038.6 1273.8 

1206.1 1195.6 2238.4 33722.3 51.58 

82 1.375 2975.8 2928.0 2550.7 2019.2 1275.3 

1209.9 1199.3 2237.7 33680.5 51.60 

83 1.625 2930.3 2884.3 2517.2 2000.3 1276.5 

1213.6 1203.8 2236.9 33640.3 51.62 

84 1.675 2803.1 2838.1 2481.8 1980.0 1277.6 

1217.1 1207.7 2236.2 33601.6 51.63 

85 1.725 2832.7 2789.1 2444.2 1958.4 1278.4 

1220.5 1211.4 2235.5 33564.5 51.65 

86 1.775 2781.7 2739.5 2406.1 1936.5 1279.0 

1223.6 1215.0 2234.8 33529.0 51.66 

87 1.825 2724.4 2683.8 2363.1 1911.6 1279.5 

1226.7 1218.4 2234.1 33495.1 51.67 

88 l.d75 2664.0 2625.1 2317.9 1885.3 1279.7 

1229.5 1221.6 2233.4 33462.9 51.69 

89 1.925 2o00.4 2563.3 2270.2 1857.5 1279. 7 

123?.. Z 1224.7 2232.7 33432.4 51.70 

90 1.975 25j2.0 2496.8 2218.9 1827.5 1279.5 

1234.7 1227.6 2232.0 33403.6 51.71 

91 2.J25 2*h2.b 2429.3 2166.7 1796.9 1279.1 

1237.0 1230.4 2231.3 33376.5 51.73 

92 2.075 2369.2 2358.0 2111.5 1764.4 1278.5 

1239.2 1233.0 2230.5 33351.1 51.74 

93 2.125 2310.4 2281.4 2052.2 1729.5 1277.7 

1241.1 1235.4 2229.8 33327.4 51.75 

94 2.175 22il.3 2204.4 1992.6 1694.3 1276.6 

1242.9 1237.6 2229.1 33305.6 51.76 

95 2.225 2147.4 2122.9 1929.4 1656.9 1275.4 

1244.5 1239.6 2228.4 33285.5 51.77 

96 2.275 2061.0 2038.9 1864.3 1618.3 1274.0 

1245. «^ 1241.5 2227.7 33267.2 51.78 

97 2.i25 1972.3 1953.2 1797.7 1578.8 1272.3 

1247.2 1243.2 2226.9 33250.7 51.79 

98 2.375 1883.5 1866.3 1730.3 1538.7 1270.4 

1248.2 1244.7 2226.2 33236.0 51.79 

99 2.425 17<;3.9 1779.2 1662.6 1498.3 1268.4 

1249.1 1246.0 2225.5 33223.1 51.80 

100 2.475 17C4.3 1692.5 1595.2 1458.1 1266.1 

1249.7 1247.1 2224.8 33212.0 51.81 

EXIT PRESSURE EXIT TEMPERATURE 

113.5 1247.7 



Pl'WEK per unit length CKW/FT) 8.04 
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1. The maximum fuel temperature is 3569.4 F and occurs at 

Z = 0.025 ft. This compares to the 3548,2 F ^at the same 
location for normal core operation. 

2. The maximum cladding temperature observed was 1279.7 F 
and occurred at Z = 1.900 ft. This compares to 1120.0 F 
at Z = 1.925 ft for normal operation. 

3. The exit sodium temperature was 1247.7 F compared to 1100.0 F 
for normal operation. 
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5.6 Listing of Code 



C 

C ' 

C VP! L^BR THikm HYDRAULICS CQDE, 

C 
C 
C 

KfAL '<FU,KCL,KBCtM,MUCIK,KU,MD(?.T,KCLl,KCL2,KFJl,KFJ2 

C 

: reN^PERATIJi^F OEPENOEMT PRfPERHES. 

KFU(T)=3. 26743-0. 2^e27F-2*Tf0.69947E-6*T**2 

KCL(T)=0.65?61EU0.91565E-2<=T-0.48934E-5*T**2+0.?6882E-8«T**3-0,16 
<793!f-ll*T**4>0.t3502E-15*T**5 

K:iC^T)-54.306-1.878E-2*'T+2.09KF-6«T*'^2 

VI$(T)= 1. 45675-0. 14807F-2*T+0.5290E-6«T»*2 

RHOi 11=59. 566-7. 95040E-3*T-0,2872E-6*1**2+0.06035E-9*T**3 

CP(T)=0. 38935-1. 10600L-4*T+?.41iar-8*T**2 

<> 

: SATURATION TEMPERATURE AS FUNCTIGM OP PRESSURE. 

C 

TSATtP) = %75541t3+0.756 7:^2*ALCC(P)-0.3'345eEl*AL0G(P)^*2K',98415«A 
^LiblP)'^=^3 

C 

C READING Ii\PUT OATA. 

C 

PEADI5»5)DFUtOCLfS,HP,H 
5 FORMATI3F10.4,F10.29FIO.I) 
R£AD(5,10)r,IN,TIN»PIN 
10 F0RMAT(E10.3f2F10.11 

PEAU(5,15)OTPOtHFtOP 
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15 F:i;MAT(El0.3,Fl0.1,I10) 

C PRINTOUT OF IMW DATA. 

C 

V*'UT^I6,20) 

I) f- JP:mT( IHlf 1X» •-''-<*»'=**'*»r£*v-^*****#t*i(i*<!*«<!!|:ijt* ' f /i3X, 

<*'i//»3X,»IflPL;T r^nti'J/] 

^i^iTE(6f25)DFUtOCU$,h,Hy 
15 FJR:iAT(3Xt'FUEL j-D. S6X, • CLAO 0.0.» ,6X, • «0D PITCHS6X, 'ACTIVF COR 
<E LcNGTHS6X,'rtIRF XP^.P PITCH 't/tiXfM fNCH )S7X,M INCH )',7X,» 
<( INCH )M2X,'( INCH )M4X,M INCH I S/,^X,F6.4,9X,F6.4,9XtF6.4, 
<HX,F6.2fl6X,F6.2,/) 

WftITE(6,30)GIN,TIN',PlN 

3> F0RM4T(3X,MNLET MASS VFLlICITY»,3X, 'CORE INLET TEHPERATJRES3X, 'CI 
<RF INLET PRESSUPF',/,5X,»( LB/HR-Fm2 IM3X,M F )M6X,M PSI 
<A )S/tTX,EI1.4,lAX,F7.2,17X,F7.2,/) 
WRITei6f35)0TPaiHE 
35 F0RHATI3X,•V0LU^'ETRIC THERMALS4Xf 'EXTPAPOLATEO HEIGHT',/, 'SOUR 
<CE STRENGTH', /,<^X,M RTU/HR-FT»*3 )MIX,MINCH ) •,/,6X,Ell.<^,HX, 
<F6.2,//) 
WRITFI6,40)0P 

'^G F0RMAT(I9X, 'OPTION', 3X, 12, ///,3X,'E NO OF INPUT DATA', 

</f3Xf '***********'^*****=ff*«*****<«5f:*** «,////) 

C 

C CONVERTING IMCHFS TO FEET. 

C 

OFU=OFU/U.O 

S=S/12,0 

0CL=DCL/I2.0 

H=H/12,0 

Ht=HE/I2.0 
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HP = HP/li,0 
PI = 3.U]59 

C 

GEOMETRY CILCULATICNS. 

r 

It 

D.v.^ = S-DCi. 

AFL=S'=*2^CnSIPI/6.0)/2,- PI/fi.0^^(0CL=^*2 + DHW^*2) 

»>.i="I-M0CL+DVJW)/2, 

Dtg=4,0*AFL/PW 

SRAT=S/nCL 

C 

C RATIO OF FRICTION FACT-i^ IN SMOOTH TRIANGULAR CORLAMT CHANNEL 

C TO CIRCULAP TUBE FRICTION FACTOR. 

C 

FRAT=-3.475i+8,0528l--S<'AT-4,7047*SPAT**2+0.91621*SftAT*«3 

C CALCULATION OF RATIO OF FRICTION FACTOR IN HIRE-WRAPPED CQOLAMT 

C CHANNEL TO THAT IN SMOOTH TRIANGULAR CHANNEL. 

C 

Y= 0,I8598-0,13640El*0Rf0,289A4El^^DF**2 
FwWM=l,195-Y«IHP-.75l/*25 

C 

C '^<ASS FLOli ^Ml CALCULATION. 

C 

Mr)!lT=GIN*AFL 

C 

C SELECTION 3F AXIAL CALCULATION INCREMENT, 

r 

f)ELL=H/100.0 

C 

C CALCULATION OF PPESSU!?r- LOSS AT INLET OF CORE. 
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U-^rClMFHOITINl/St-CO.O 

f'E=Ue*DEQ=^PHG(TIM)/VIS(Tirj)*3600.G 

FCIPC=0.1fl4/RE**C.2 

LiELPF-F*=:j£LL*;<Ht]( TIN )^UB**2/nF0/2.. 732. 2*23.0 

OELP-=0.5*RHOlTIN)*U3**2/2,0/32.2 

WRITE(6,'f2) 

FnRMAT(5X,MT THE INLFT OF THE CCPF',//) 

lv'^ITE(6,4i))UB,RF,FCIP.C 

t.*j F0RMAT(3X, 'CHANNEL AVE. VELOCITY S 2X, 'REYNnLOS #S3X, 'FRICTHM FA 

<CT0R',/,7X,MFT/SEC)',/»3X,F6,2,12X,F8.l,8X,F7.4,//l 
KRITE(6,50I 

c . 

C CALCULATION OF CONDITIONS AT BEGINNING OF FIRST INCREMENT. 

C 

Pl=PIN-IOELPFtDELPE)/lH.O 

T3l=TIN 

1 = 1 

> C CALCULATIQNJ OF AXIAL CCOPCINATF OF CENTER OF FIRST INCREMENT. 

f 

ZC=-H/2.0+nELL/2.0 

C 

C ??.INT OUT OF COLUMN HEADRiGS. 

C 

WRITEIo,55) 

•55 F0RMAT(2Xt*NnS4X,»7-Lrr.S6X,»T0',7X,'Tl/4SSX,»Tl/2'i6X,»T3/4S7 
<X,*TFS6X,'T CLAD't4X,'T :IULKS4X,»T SAT«,3X,'H T COEf- ' ,4Xt » OELTA 
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<PS//) 

c 

C CALCULATIUN CF AVEi^AGE BULK TEMPERATURE FOR INCREMENT.' 

r 

\^ 

60 OTPAV=QTPC*COS(PI*ZC/HE) 

T!^2=TPl+0TPAV*PI*DFU**2.0*DELL/4,0/M00T/CP|TBU/2.0 

I 

C THE DIVISHN (^Y ? IN THH PREVIOUS CALCULATION IS REQUIRED BF. - 

C CAUSE ONLY HALF OF THE ENERGY GEMERATEO IN THE FUEL ROD ENTERS 

C THE CHANNEL, 

C 

TBAV=(TFI+TP2)/2.0 

C 

C CALCULATION OF PRESSURE DROP A^'D AVERAGE PRESSURE FOR INCREflENT 

C 

tirGIN/PHQ(TaAV)/3600.C 

l-^J3«l}E0«RHa(TyAV)/VIS(T3AV)*3600.0 
FCIRC=0.184/RE**0.2 
F=FCIRC*FRAT*FWWM 

OtLPF=F*DELL»RH0(TRAV)*ll3**2,0/DEO/2. 0/32,2 

P2=Pl-mPF/144.G 

PAV=(Pr+P2l/2,0*144,0 

C 

C CALCULATION! OF AVERAGE SATURATION TEMPERATURE FOR INCREMENT. 

C 

TSATAV=TSAT(PAV) 

C 

C CALCULATIOM OF HEAT TRANSFER COEFFICIENT FOR INCREMENT. 

C 

PR=CPiTflAV)*VIS(T8AV)/KRC{TeAV) 
PE=RE*PR 

NU=-88,0+U0.e*SRAT-4U6r=^SPAT**2 
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PPG=3'36.J^^SRAT-110.0 
IF(PF-Pt-n)70,70,t5 

70 HTC=^iU-^Kt3C(T3AV)/OEO 

r 

C CALCULuTIO'^ uF TFMPERATURF AT UUTSID^ SURFACE OF CLADDING FOP 

C L^JC^X.-ltM. 

C 

TSi\V=TB'\V+0TDAV*0FU'^*2.0/4.0/Dr,L/HTC 

C 

C CHECK F0(^ PQSSieiLTTY GF SURFACE BCILING. IF DETECTEDi "SUR- 

C FAC; BPILIMG" PRINTED OUT AND CALCULATIONS TERMHATED. 

C 

IF(TSAV-TSATAV)95,75,75 
75 WRITF(^,f)Cl 

8'J FORMAT!/, lOX,' SMPFACE BOILING',/) 
WRIT5(6,85) 

FORMAT!/, lOX, 'SURFACE TEMPERATURE », 1JX,» SATURATION TrMPERATJRE' ,/) 
W<^ITE(6,90)TSAV,TSATAV 
90 I- jRMATll6X,F9,2,lOX,F9.2) 

00 Tn 160 ;i . . 

C . 

C CALCULATION^ OF Tc.MPfRATUfE AT rUTSIDE SURFACE OF FUEL FO^ IN- 

C CKtMENT. DETERMINED BY ADDING TEMPERATURE RISF ACROSS CLAD0IN5 

C TO TFMPFPATURE AT OUTSIDE SURFACE QF CLADDING. AN ITERATION 

C TECHNIQUE IS USED TC OBTAIN THF CLADDING THERMAL CONDUCTIVITY 

C TO WITHIN I?. 

C 

9!) KCL1=KCL{TSAVI 

lOU TF=TSAV+QTPAV<'DFU**2.0/8,0/KCLl*ALnG(OCL/DFGI 
TCLAV=(TSAV+TFI/2.0 
KCL2=KCL(TCLAVI 
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FRft=AeSI{KCLl-KCL2)/KCL2) 
Ifit^P-O.JODUCllOTlOS 
lOi) KCLUKCL2 
G;j TO 100 

C CALCULATION ff TE^PFRATiJPE AT CnMTFR OF FUEL FOR INCREMENT. AM 

C ITHATIGM TECHNIOUC IS USED TO •'^PJAIN FUEL THERMAL CONDUCTIVITY 

C ^ITHIN 1%. 

ah 

t 

UO KFU1=KFUI^F) 

115 TG=TF+QTPAV*0FU^*2.0/U.0/KFU1 

TFAV=(TF+TC)/2.0 

KFU2=KFU(TFAVI 

ERR=A6S{(KFUl"KFU2)/2,0) 

IFIEKR-O. 001)125, 125tl20 
120 KFU1=KFU2 

GO TO 115 

C 

C CALCULATIUN OF FUEL TE'^PERATURE AT 1/4, 1/2, AND 3/4 FUEL RADII 

C m INCJ^EMENT, 

C 

125 T14=TO-QTPAV*DFU«*2.2/16.0/KFU1/16,0 
T12=T0-gTPAV*DFU**2,0/l6.0/KFUl/4.0 
T54=T0-0TPAV*0FU**2.0/16.0/KFU1/16,0*9,0 

C 

C f^HINTOUT OP RESULT FOR CALCULATION INCREMENT, 

C 

WRITcI6,130}I,ZC,T0,T14,T12,T34,TFtTSAV,TBAV,TSATAV,HTC,0ELPF 
no FORMAT(1X,I3,F10.3,9F1C.1,4X,F6.2} 

c 

C CHECK TO SEE IF CALCULATIONS HAVE BEEN COMPLETFD FOR EMTIRE 

C CORE. 
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f 

IF(I-ljC)l35,U0,U0 

C 

C iETTlHG CllN'OITIONS MIR tlEGINMNG OF NEXT CALCULATION INCPHMENT. 

C 

U^i 1 = 1 + 1 

K=7X+DELL 

T81=TB2 

Pl=P2 

r 

C TPANSFFPRUG TO START CALCULATIONS FOR NEXT INCREMENT. 

C 

GO TO 60 

C 

C CALCULATION OF CORE EXIT CONDITIONS OF COOLANT. 

C 

HC UBsGIK/PH0ITB2l/3600.0 

R6=UB*0£Q*RH0(TB2)/VIS(TB2)*3600.0 

FCIRC=0.I84/PF»*0.2 

F=FCIRC*FRAT*FWWM 

0ELPF=F*DELL*RH0(TB2)*UB**2/DE0/2. 0/32. 2*103.0 
D£LPE=RH0{TB2)*U6**2/2.0/32.2 

C 

C CALCULATION OF CORE FXIT CONOITIONS OF COOLANTS . 

r 

PfX=P2-(OELPF+0eLPFI/U4.0 
T£X=TB2 

r 

C PRINTOUT OF CORE EXIT COOLANT CONDITIONS. 

C 

wRITE:(6,l45) 

145 FOKMAT(/tICX,'FXIT PPESSURl:MOXt»EXlT TEMPERATUPF' ,/ ) 
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WKlTEl6,150)PrX,T^X 
15) f j;i^'U(UX,P6.ltl8X,F6,l) 

C 

C rM.a'UTirj\ OF AVEi^AG'': Km or FUEL m IN KW/FT. 

C PMis'TJl'T OF AVe^AGF ^^'■i'r.K 

r 

i-.'VL=OTPr/-Hf.^r.Fll*rt2/2.0:^SlM(H=^h/2.0/HE)/H 
>KlTE(6,l5^.)PrVL 

155 F3R:^AT(//,5X,' PGwER PE^ JNIT LErGTH (K'/^/FT) ' ,Fl?. 2, /////) 
160 STQo 

I 

00 

I 
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7.0 Problems 



1. By performing an energy balance on a differential element with sides dx, 
dy, and dz, derive the steady state heat conduction equation in Cartesian 
coordinates • 

2. Derive the steady state heat conduction equation in cylindrical coordinates. 

3. Develop an expression for the temperature distribution in a long, thin, 
hollow cylindrical fuel element with inner radius r^, outer radius r^, and 
volumetric thermal source strength q"'. Assume that the inner surface is 
insulated and at temperature T^. 

A. A long, thin, hollow cylindrical fuel element with volumetric thermal source 
strength q'" is cooled on both the inner and outer surfaces (r « r^ and 
r = r surfaces) such that the surface temperatures are and T^, respective- 
ly, fierive an expression for the location of the maximum^ radial temperature 
in the fuel. 

5. Calculate the ratio of peak power to average power for a fuel rod of constant 
cross section assuming that the volumetric thermal source strength varies as 
q'" = q cos (nz/Hg) and that the extrapolation length is 10% of the core 
height. ^Compare this ratio to that for the case where the extrapolation 
length is assiamed to be zero. 

6. Show that for a cylindrical shell without internal heat generation the use 

' of an average constant thermal conductivity evaluated at the arithmetic mean 
of the inner and outer surface temperatures is equivalent to a linear 
variation of thermal conductivity with temperature. 

7. A LMFBR fuel rod consists of 0.250 in. diameter UO2-PUO2 fuel pellets in a 
0.250 in. ID, 0.300 in. OD tube of stainless steel cladding. The heat trans- 
fer coefficient at the ou* ^.r surface of the fuel rod is 30,000 Btu/hr-ft -F. 
The thermal conductivities of the fuel and cladding are 1.95 Btu/hr-ft-F and 
12.8 Btu/hr-ft-F, respectively. The coolant bulk temperature is 900 F and 
the volumetric thermal source strength is 1.8 x 10^ Btu/hr-ft^. Determine 
the maximum fuel temperature and the maximum cladding temperature. 

8. Same as Problem 7 except the ID of the cladding is 0.270 in. The 0.020 in. 
diametral gap between the fuel and cladding is filled with helium gas. 
Assuming that the heat transfer through the helium is by molecular conduction 
only (k, - 0.140 Btu/hr-ft-F) determine the maximum fuel and cladding 
temperatures. 

9. A four foot long LMFBR fuel assenbly contains 0.260 in. CD fuel rods spaced 
0.320 in. in a triangular array. The wire wrap helical pitch is 9 inches. 
The average conditions in the core are 150 psia and 950 F. The mass velocity 
of the liquid sodium is 5.0 x I06 Ibm/hr-ft^. Determine the pressure drop 
through the core. 
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10. Determine the effect on pressure drop of a 10% increase in cladding CD 
for the fuel in Problem 9. Assume the spacing, coolant conditions and 
mass velocity remain the same. (Note the wire wrap diameter must be 
reduced by 0.026 in. to maintain the same spacing of the fuel rods). 

11. Determine the heat transfer coefficient for the fuel and flow conditions 
of Problem 9. 

12. A LMFBR fuel consists of 0.250 in. diameter UO2-PUO2 fuel pellets contained 
in 0.250 in. ID, 0.285 in. CD type 347 stainless steel cladding. The fuel 
rods are spaced 0.350 in. in a triangular array. At a particular location 
in the core the pressure is 100 psia, the liquid sodium temperature is 
1000 F and the coolant velocity is 24 ft/sec. The volumetric thermal 
source strength is 1 x IO8 Btu/hr-ft^. Determine the heat transfer 
coefficient. 

13. Determine the coolant velocity required to keep the temperature rise across 
the convective layer below 12 F for the fuel and coolant conditions of 
Problem 12. 

14. In a nuclear fuel rod the volumetric thermal source strength varies axially 
as q'" (z) = ^o' " (IIZ/H ). The energy generated in this rod is 
removed by coolant channels w^ich have a combined mass flow rate of m and 
an inlet temperature of T^^ Designate the fuel diameter as D^ and the 
cladding CD as D . Develop'an expression for the axial variation of the 
temperature at tfie inside surface of the^ cladding. Your expression should 
contain only T^^, q^'", D^, D^, H, H^, m, C^, h, k^, and Z. 

15. From the expression developed in Problem 14 determine the location of the 
maximum cladding temperature in the reactor. 

16. Determine the location of the maximxxra fuel temperature for the fuel rod- 
coolant channel system of Problem 14. 

17. A 1400 Mw(t) liquid sodium cooled LMFBR contains 25,000 fuel rods each 
composed of 0.200 in. diameter UO2-PUO2 fuel pellets in 0.230 in CD tubes 
of type 347 stainless steel cladding. The fuel rods are spaced 0.290 in. 
in a triangular array. The reactor core height is 60 in. The coolant 
inlet conditions are 250 psia, 700 F and the average mass velocity is 5.4 
X 10^ lbm/hr-ft2. The thermal source strength varies as q"' = ^o' ' ' 
(IlZ/H ) and H = 65 in. Determine the maximum fuel temperature and the 
maximum claddfng temperature in the core. 

18. Determine the % reduction in coolant flow that will cause the maximum cladding 
temperature to approach the allowable limit of 1250 F for the reactor described 
in Problem 17. 

19. Determine the % overpower that the reactor described in Problem 17 can produce 
without exceeding the 1250 F maximum cladding temperature limit. What is the 
coolant temperature at the core exit for this limiting case? 
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20. A 1200 Mw(t) liquid sodium cooled LMFBR contains 30,000 fuel rods with 
an active fuel length of 53 in. Each fuel rod consists of 0.180 in. 
diameter U02'PviO2 pellets enclosed in 0*210 in. OD type 347 stainless steel 
claddingr The fuel rods are spaced 0.260 in. in a triangular array. The 
inlet cooJ.ant conditions are 700 F, 250 psia and the average velocity in 
the core is 32 ft/sec. The axial power distribution is a truncated cosine 
function. Assume an extrapolation length of 3 in. Determine the location 
and magnitude of the maximum fuel, cladding, and coolant temperature. 

21. Determine the % of overpower that would result in a 1200 F maximtim allowable 
cladding temperature limit being exceeded for the core described in Problem 
20. At what axial location does this 1200 F temperature occur? 



93 



